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Abstract 


The  contents  of  this  report  covers:  (i)  development  of  optimal  geometry 
for  crowned  spur  gears;  (ii)  methods  for  their  generation;  and  (iii)  tooth 
contact  analysis  (TCA)  computer  programs  for  the  analysis  of  meshing  and 
bearing  contact  of  the  crowned  spur  gears. 

The  developed  method  for  synthesis  is  used  for  the  determination  of  the 
optimal  geometry  for  crowned  pinion  surface  and  is  directed  to  reduce  the 
sensitivity  of  the  gears  to  misalignment,  localize  the  bearing  contact,  and 
guarantee  the  favorable  shape  and  low  level  of  the  transmission  errors. 

A  new  method  for  the  generation  of  the  crowned  pinion  surface  has  been 
proposed.  This  method  is  based  on  application  of  the  tool  with  a  surface  of 
revolution  that  slightly  deviates  from  a  regular  cone  surface.  The  tool  can 
be  used  as  a  grinding  wheel  or  as  a  shaver.  The  crowned  pinion  surface  can  be 
also  generated  by  a  generating  plane  whose  motion  is  provided  by  an  automatic 
grinding  machine  controlled  by  a  computer. 


The  TCA  program  simulates  the  meshing  and  bearing  contact  of  the 
misaligned  gears.  The  transmission  errors  are  also  determined. 
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Spur  years  are  widely  used  to  transmit  mechanical  power  between 
parallel  axes  of  drive  shafts.  Namely,  such  gears  are  used  in 
planetary  trains  of  helicopters.  It  is  well  known  that  spur  involute 
years  are  very  sensitive  to  gear  misalignments.  Misalignment  will 
cause:  the  shift  of  the  bearing  contact  toward  the  edge  of  the  gear 

tootn  surfaces  and  transmission  errors  that  increase  gear  noise.  The 
previous  metnods  for  crowning  of  pinion  tooth  surface  have  been 
directed  at  providing  of  a  favorable  bearing  contact.  The 

transmiss ion  errors  of  the  misaligned  years  have  been  ignored.  The 
new  approach  developed  oy  the  authors  of  this  report  is  directed  at 
the  syntnesis  of  spur  involute  gears  with  a  crowned  pinion  surface 
taut  provides:  (1)  a  parabolic  type  of  low  transmission  errors  for 
tne  misaligned  years  and  (ii)  a  localized  bearing  contact.  New 
metnods  for  the  generation  of  crowned  pinion  surface  proposed  by  the 
autnors  can  be  readily  applied  by  the  gear  industry. 
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1.  Introduction 

Spur  involute  gears  are  very  sensitive  to  year  misalignments  that 
cause:  (i)  the  shift  of  the  bearing  contact  to  tne  edge  of  tne  gear 
tooth  surfaces  and  (ii)  transmission  errors  that  increase  gear 
noise.  Many  efforts  have  been  done  to  improve  the  bearing  contact  of 
misaligned  spur  gears  by  crowning  of  pinion  tooth  surface.  Wildhaber 
(Ivor)  nas  proposed  various  methods  of  crowning  that  can  be  achieved 
in  the  process  of  year  generation.  Maay  engineers  (see  References) 
have  usee  erownin_,  as  (1)  longitudinal  corrections  (tig.  1.1a),  (ii) 
mouirieu  involute  tooth  profile  uniform  across  facewidtn  (Fig.  l.lo), 
(in;  combination  of  longitudinal  correction  and  uniform  modified 


prolife  (tig.  l.lc)  and  (iv)  topological  mod i f ica t ion  (Fig.  l.ld)  that 
can  proviso  any  deviation  of  the  crowned  tooth  surface  from  a  regular 
involute  surface. 

ihe  main  purpose  of  the  previously  proposed  methods  for  crowning 
was  to  improve  the  bearing  contact  of  the  misaligned  years.  This 
resulteo  in  only  part  of  tne  solution.  The  transmission  errors  of  the 
misaligned  spur  gears  (a  main  source  of  the  noise),  were  ignored.  The 
influence  01  gear  misalignment  on  the  transmission  errors  has  not  been 
investigated .  Maay ' s  method  of  topological  mod i f i ° at  ion  does  not 
quantity  the  relationship  between  the  surface  deviation  and  tne 
t runsmiss  ion  error.  Also  the  optimal  geometry  for  the  crowned  pinion 
s  ;  r Lace  nas  not  Deen  proposed. 

ihe  contents  ul  this  report  cover  tne  solutions  to  the  following 
t,  r  j.j  ie..,s :  (1)  optimal  geometry  of  a  crowned  pinion  tooth  surface; 
lii)  now  metnou  or  irowniig  based  on  application  of  a  tool  with  a 
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tooth  surface  is  based  on  the  following  considerations: 

(i)  Misaligned  spur  gears  with  a  crowned  ninion  tooth  surface 
can  provide  transmission  errors  A  2  ( *  ^ )  of  two  types  that  are  shown  in 
Fij.  1.2a  and.  1.2b,  respectively.  The  transmission  errors  are 
ieterr.ined  with  the  eouation  as 


: (  :  1 }  =  *2{*i: 


%  *1 


(1.1) 
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ear  rotation; 


*  2  (  y  )  is  the  function  that  relates  the 
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gears  are  misaligned?  =  —  $  i  is  the  theoretical  relation  between 

the  angles  of  rotation  of  the  gears  in  the  ideal  case  where  the  gears 
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errors  of  type  1  (Fig.  1.2a).  This  is  the  reason  why  the  authors  have 
decided  to  synthesize  a  specific  crowned  pinion  tooth  surface.  Such  a 
pinion  can  provide  transformation  of  rotation  with  a  parabolic  type  of 
transmission  error  function  when  it  is  in  meshing  with  a  gear  with  a 
regular  involute  tooth  surface.  This  type  of  function  of  transmission 
errors  is  synthesized  for  ideal  years  that  do  not  have  any 
misalignment.  Then,  the  tendency  to  provide  parabolic  transmission 
errors  can  be  extended  for  the  misaligned  gears  and  the  discontinuance 
on  mesnin.,  can  be  avoided.  it  is  necessary  to  empnasize  that  the 
proposed  method  of  synthesis  provides  a  shape  in  the  middle  cross- 
section  of  tne  pinion  tootn  that  deviates  in  a  certain  way  from  the 
involute  curve  that  is  shown  in  Fig.  1.3.  The  longitudinal  deviation 
trom  a  straight  line  is  not  related  witn  the  transmission  errors  but 
with  the  desired  dimensions  of  the  instantaneous  contact  ellipse  for 
tne  jcar  tootn  surfaces.  The  proposed  pinion  tooth  surface  can  be 
generated  by  a  plane  as  the  chosen  generating  surface.  The  motions  of 
tne  plane  with  the  respect  to  pinion  must  be  controlled  by  a  computer. 

(lii)  Another  method  of  pinion  crowning  is  based  on  application 
ol  a  surface  of  revolution  that  slightly  deviates  from  a  regular  tool 
conical  surface  (Fig.  1.4).  Such  a  tool  can  be  used  as  a  grinding 
Wiicei  or  as  a  saver.  i’iie  i, lotions  of  the  tool  and  the  gear  being 
generated  are  related  similarly  to  the  motions  of  a  rack-cutter  and 
Luc  jC-ur  (see  Chapter  4).  A  tool  with  a  regular  conical  surface  can 
jene  rate  a  crowned  pinion  tooth  surface  whose  middle  cross-section 
represents  an  involute  curve  (Fig.  l.J).  However,  this  type  of  a 
crowned  pinion  tooth  surface  is  not  desirable  because  the  misaligned 
jeers  r,r ovide  transmission  errors  of  type  1  (Fig.  1.2a).  This  is  the 


J 


reason  why  a  cool  with  the  surface  of  revolution  is  to  be  used  instead 
of  a  conical  surface. 

(iv)  The  evaluation  of  the  bearing  contact  and  transmission 
errors  for  the  misaligned  gears  as  well  as  the  investigation  of  the 
influence  of  errors  of  gear  assembly  requires  application  of  TCA 
programs.  Such  programs  have  been  developed  by  the  authors.  The 
programs  are  eased  on  the  following  algorithms: 

id)  The  contacting  gear  tootn  surfaces  are  represented  in  a 
iixed  coordinate  system,  Sc,  tnat  is  rigidly  connected  to  the  gear 
nous  in_,  ( c  ig  .  I .  b  )  . 

(u)  Tnc  continuous  tanyencv  of  gear  tooth  surfaces  is  provided 
11  tne  pos i t ion-vectors  and  surface  unit  normals  for  the  contacting 
sartaces  coincide  ac  the  contact  point  at  any  instant.  Then,  we  are 
uoie  to  determine  the  path  of  contact  on  the  gear  tooth  surfaces  and 
mo  relations  between  tne  angles  of  rotation  of  the  output  and  input 
gears.  Knowing  function  ^  (4^),  we  can  determine  the  deviations  of 
tnis  function  from  the  prescribed  linear  function  i.e.  the 
transmission  errors. 

(c)  bue  co  tne  elasticity  of  the  gear  tooth  surfaces  the  surface 
contact  is  spread  over  an  elliptical  area.  The  dimensions  and 
orientation  of  the  instantaneous  contact  ellipse  depend  on  the 
principal  curvatures  and  principal  directions  of  the  contacting  tooth 
suriaces.  The  ocarin,,  contact  is  determined  by  the  developed  TCA 
program  as  tne  set  of  the  contact  ellipses  that  move  over  the 
contact  in j  surfaces  in  tne  process  oi  motion. 


2.  Optimal  Pinion  Geometry 

2.1  Basic  Consideration 

It  is  assumed  that  only  the  pinion  will  be  provided  with  a 

crowned  tooth  surface  while  the  year  is  provided  with  a  regular 
involute  surface.  The  crowned  pinion  tooth  surface  deviates  from  a 
regular  involute  surface.  The  topoloyy  of  the  crowned  pinion  tooth 
surrace  must  satisfy  the  following  requirements:  (i)  the  pinion  tootn 
surface  is  a  continuous  one;  (ii)  the  pinion  and  year  tooth  surfaces 
are  in  contact  at  a  point  at  any  instant;  (iii)  the  principal 
directions  and  curvatures  of  the  contacting  surfaces  must  be  related 
in  sucn  a  way  that  tne  instantaneous  contact  ellipse  will  be  of 

appropriate  dimensions;  (iv)  tne  function  of  transmission  errors  of 
tne  misaligned  years  snould  be  of  a  parabolic  type  and  its  level  must 
not  extend  the  prescribed  limits.  The  satisf ication  to  the  described 
conditions  is  considered  in  the  following  sections  of  this  part. 

2.2  Parabolic  Type  of  Function  of  Kinematic  Errors 

Initially  we  will  consider  that  the  gears  are  not  misaligned,  the 
pinion  tooth  surface  is  crowned  and  the  gear  tooth  surface  is  a 

regular  involute  surface.  The  tooth  surfaces  are  in  point  contact  at 
every  instant  and  tne  points  of  contact  lie  in  the  middle  cross- 
section.  Thus,  initially  we  consider  tne  meshing  of  year  tooth 
surfaces  as  a  planar  gearing.  however,  when  simulating  the  meshing  of 
gears  og  itrt  program,  we  will  consider  tne  mesh  of  them  as  spatial 

j  e  a  r  s  . 

Figure  s.z.L  snows  the  shapes  of  the  gears  which  are  in  contact 
at  a  point  that  lies  on  the  center  distance.  The  shape  of  the  gear 


.v.v„vv)v.v 


. Hi  « 


tooth  profile  is  a  regular  involute  curve,  and  the  shape  of  the  pinion 
tooth  profile  is  to  be  determined.  It  is  given  that  the  shapes  must 
provide  transformation  of  rotation  represented  by  the  function  (Fig. 
2.2.2a) 


> 2  ^  $ i )  4*]_  A<j>2  (  -^ ) 


(2.2.1) 


where  A o  ^  *»* i )  i-s  a  function  of  a  parabolic  type  that  represents  the 

function  of  transmission  errors;  and  N2  are  the  numbers  of  gear 

teeth.  We  have  to  emphasize  that  the  proposed  method  for  synthesis 

provides  low  transmission  errors  for  aligned  gears  and  the  function  of 

transmission  errors  is  a  parabolic  one.  Then,  the  misaligned  gears 

will  keep  the  tendency  to  transform  rotation  with  the  same  type  of 

transmission  errors  and  the  discontinuance  of  contact  by  meshing  will 

be  avoided.  This  means  that  we  have  to  avoid  the  appearance  of 

transmission  errors  that  are  shown  in  Fig.  1.2a.  Function 

N1 

4, 2  ( 4>  i )  contains  a  linear  function  represented  by  jj-  <j>^  that  relates 
the  angles  of  rotation  of  the  "ideal"  input  and  output  gears.  Figure 
2.2.2a  shows  the  function  2  (  <f>  ^ )  •  The  function  of  transmission 
errors  is  represented  to  a  larger  scale  in  Fig.  2.2.2b  as  a  periodic 
function  of  a  parabolic  type.  The  derivative 


A  9  (  ^  ^  )  )  ^  ^  (  A<J>  2  ^  4*  1  ^  ) 


(2.2.2) 


is  shown  in  Fig.  2.2.2c.  For  a  "pure"  parabolic  function  this 
derivative  is  a  linear  function.  It  is  evident  that  the  second 
derivative  represented  by 


— 2  (  A  4>  ^  (  <f>  ^  )  )  -  - 2 

d  (J)^  d  (J)-^ 

is  negative.  The  transformation  of  rotation  by  the  gears  will  be 
performed  with  a  non-constant  gear  ratio 


m 


21 


"2 

“1 


A  <J)  2  (  ^  ^  ) 


(2.2.3) 


The  instantaneous  center  of  rotation  I  is  moved  along  the  gear  center 
distance  in  the  process  of  gear  rotation.  Figure  2.2.1  shows  the 
instantaneous  location  of  I. 

More  details  about  tne  yeometry  of  the  pinion  that  corresponds  to 
tne  discussed  function  <f>^(<f>^)  will  be  given  in  Chapter  3. 


z . 3  Relation  Between  the  Principal  Curvatures  and  Directions 

The  relation  between  the  principal  curvatures  and  directions  of 
tne  contacting  surfaces  and  the  derivative  of  the  gear  ratio  is 
represented  as  follows  (Litvin,  1987) 


b 

li 

D  1 2 

b13 

b 

21 

b22 

b23 

=  0 

b 

31 

b32 

b  33 

For  tne 

case 

ot 

spur  gears  with 

the 

o  leiaen  ts 

of 

the 

determinant  of 

Eq . 

tol lowing  equations 


(2.3.1) 

crowned  pinion  tooth  surface  the 
(2.3.1)  are  represented  by  the 
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U  2 1  D 1 2 


bll  b12 


b 1 2  b  22 


b  3 1  b  1 3 


-KU)  +  K(2) 

I  I 


(1)  .  (2) 
KII  +  KII 


(2.3.2) 


(2.3.3) 


°  1 3 

ii 

b0  . 

2  j  -J 

-e  j  •  e 


(2.3.4) 


+  (  2  1 5  )  2ra  (  nxk2  ) .  (rj^-C) 


(2.3.5) 


Figure  2.3.1  shows  the  tangent  plane  to  the  gear  tooth  surfaces  at  the 
point  of  contact  I  (see  Fig.  2.2.1).  The  principal  directions  are 
represented  by  the  unit  vectors  e^  and  e^  that  are  the  same  for  the 
uotn  contacting  surfaces.  Vector  eT  is  shown  in  Fig.  2.3.1  and  is 
represented  in  coordinate  system  3f  as  follows 


sing 

c 

-COS  g 


( 2.3.6) 


wiiere  is  tue  pressure  an^le  at  point  I.  Vector  is  directed 

parallel  to  the  negative  z^-axis  and  is  represented  by 


u 

-1 


(2.3.7) 


The  surface  unit  normal  vector  is  represented  by  (Fig.  2.2.1) 


~f  ~IX~II 


COS  Ij) 

sin^ 


(2.3.3) 


Tne  angular  velocities  of  rotation  are  represented  by  (Fig.  2.2.1) 


(  1)  _  ( 1) 

(0  -  10 


(2)  _  (2) 
0)  -  to 


(2.3.9) 


Vectors  and  represent  the  velocities  of  the  contact  point  in 

transfer  motion,  while  gears  rotate.  Thus  (Fig.  2.2.1) 


(1)  n  j 

to  I 


(2)  _  (2)  n 
~tr  "  “  °2 1  0 


(2.3.10) 


Vector  C  and  r^  are 


V  -  1 


r(  =  01I  1 


(2.3.11) 


mm 


where  0,0,  =  C°  =  +  ~  ,  0,1  =  ~  0oI  =  — ;  p  is  diametral 

12  zp  2p  1  2p  2  2p 

pitch  and  C°  is  ideal  or  nominal  central  distance. 


0,1  =  ^ 
1  2p 


< j  and  are  the  first  and  the  second  principal  curvatures  of  the 

pear  tooth  surface  represented  by 


.(  2)  _ 


0  -,1s  m  a 
2  rc 


(2)  _  n 
KII  "  0 


(2.3.12) 


wnere  u  ,1  =  is  the  radius  of  tne  gear  pitch  circle.  The  positive 
sign  of  tne  curvature  indicates  that  the  curvature  center  lies  on  tne 


positive  direction  of  the  surface  unit  normal.  The  unit  vector  k,,  in 


r-(i «  (2.3. 5)  coincides  with  the  unit  vector  of  w  .  Thus  k.->  =  where 
k  is  tne  unit  vector  of  the  z^-axis  (Fig.  2.2.1). 

Equations  (2.3.1)  -  (2.3.12)  yield  the  following  relation  between 
tne  derivative  m'^  =  (m2i($i))  and  the  principal  curvatures  of  the 
crowned  pinion  tooth  surface 


1 4  ^  tan 


_ ^P 

h  x  5  i  n  1(4  c 


2p  (  lVK  I  -  NlSin*c) 
N  “s  1  n  +  VN7“‘ 


(2.3.13) 


Inc  jOct  1  is  to  provide  a  paracolic  type  of  the  transmission  error 


function  at  least  in  the  neighborhood  of  the  main  contact  point  I 


This  requirement  can  oe  satisfied  with  the  negative  value  of  mo¬ 
using  Eg.  (2.3.13),  we  can  control  the  value  of  (it  is  negative) 

oy  varying  the  values  of  It  is  evident  that  is  negative  if 

.  Here:  rr—fi— —  is  the  curvature  at  point  I  if  the 
pinion  would  oe  provided  wich  a  regular  involute  curve;  <1^  is  the 


Using 

Eq. 

(2.3. 

oy  varying 

the 

K(1) 

|  > 

2p 

N  .  s  i 

pinion 

V/Ould  D€ 

curvature 

of  t. 

cross-section.  We  have  emphasized  that  the  pinion  tooth  surface  must 
oe  deviated  in  the  longitudinal  direction  (Fig.  1.3)  to  provide  the 
point  contact  of  tne  mating  tooth  surfaces  and  localize  the  bearing 


contact . 
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3.  Design  of  the  Crowned  Pinion  Tooth  Surface  for  the 
Prescribed  Transmission  Error  Function 

3.1  Basic  Concept 

Consider  that  the  function  & 2 ( 4, ^ )  that  relates  the  angles  of 
rotation  of  the  gears  is  represented  by  (Fig.  3.3.1). 


*2ui}  “  +1  jq  +  brai+i 


(3.1.1) 


Hero  : 


—  is  the  linear  part  of  function 


4>  2  ( 4>  1 )  and  A<>2  (<!]_)  =  D  i~  ax<l>4  is  prescribed  parabolic  function  of 
transmission  errors. 

The  gear  is  provided  with  a  regular  involute  surface.  If  the 
gears  are  aligned,  the  crowned  pinion  tooth  surface  and  the  gear  tooth 
surface  are  in  mesh  in  the  middle  cross-section.  The  contact  of  the 
year  tooth  surface  can  be  considered  as  the  contact  of  two  planar 
shapes.  Our  goal  is  to  determine  the  shape  of  the  pinion  tooth 
surface  in  the  middle  cross-section. 

Ue  consider  that  the  transmission  errors  in  the  region  of 
-  ~  s  are  on  average  equal  to  zero.  Thus 


[  b  .  —a  ..  <jj ,  )  d  (j, .  =  0 


(3.1.2) 


equation  (3.1.2)  yields  that 


(^)2 

d\ 


(3.1.3 


Henceforth  we  will  use  the  following  coordinate  systems 


(Fig.  3.1.2):  sl(xl,Yl)  and  S2(X2,Y-^  that  are  rigidly  connected  to 

the  pinion  and  gear  respectively;  and  Sf(Xf,yf)  that  is  the  fixed 


coordinate  system  employed  in  the  housing, 


The  determination  of  the  pinion  tooth  shape  is  based  on  the 


following  considerations: 


(i)  The  gear  tooth  shape  e2  is  a  regular  involute  curve 
represented  in  coordinate  system  S2  by  the  vector-function  ^(Aq); 
the  shape  unit  normal  is  represented  by  n 2  (  4>q  )  • 


(ii)  using  the  coordinate  transformation  in  transition  from  S2 


to  Sf,  we  can  represent  the  family  of  curves  in  Sf  as  follows 


[r^2)]  =  [Mf2: 


^  ^Lf2-*  I‘n21 


(3.1.4) 


(iii)  The  equation  of  meshing  is  represented  by 


nJ2).V^21)  =  0 

~f  ~f 


(3.1.5) 


(21) 

where  is  the  sliding  velocity  at  the  contact  point. 


bquations  (3.1.4)  and  (3.1.5)  yield  the  following  relation 


f  ( c i  G  f  A  2 )  -  0 


(3.1.6) 


that  is  called  the  equation  of  meshing. 


(iv)  The  following  equation 


r  g  ]  [  r  ^  2  ^  ] 
1  If J  1  f  J 


lMif ]  2)  [ r  2  ) 


(3.1.7) 


•  • 


and  the  equation  of  meshing  (3.1.6)  determine  the  shape  of  the  pinion 


tooth  in  tne  middle  cross-section.  The  surface  of  the  pinion  can  be 
determined  based  on  the  shape  in  its  middle  cross-section. 

3.2  Derivation  of  Pinion  Tooth  Surface 
tear  Tootn  Involute  Shape  Equations 

Tne  ^ear  involute  shage  and  its  unit  normal  are  represented  in  S2 
as  toilows: 


=  r2  lsin$G-$Gcosij>ccos(  *G-0C)  J 

y.  =  r.  (c(js$r+<f)  ,cosi(i  sin(  <pr~ip  )  ]  (3.2.1) 

Z  Z.  w  O  O  v3 


z  ,  =  0 


n 


2 


-cos  (  ) 

O  c 

sin ( ) 
0 


(3.2.2) 


He  re  : 

rz  is 

tne 

raoius  of  the 

gear 

an.jle 

cillG  'P  . 

is 

the  involute 

curve 

tUOtil 

S  n  cLjJ  G  • 

ilatrix  p-lf 2 J  is  represented  oy 


pitch  circle,  \pc  is  the  pressure 
parameter;  n2  is  the  normal  of 


l  i-i 


f  ;  J 


-  c  o  s  _  s  i  n  $  U  0 

o 

-sin<j>„  -cos<j>0  U  L 


(3.2.3) 


Us  my  Eqs .  (3.1.4),  (3.2.1),  (3.2.2)  and  (3.2.3),  we  represent  the 

family  of  shapes  £ ^  and  the  shape  unit  normal  in  Sf  as  follows 


■  (  2)  , 
r£  1  "  c2 


-sin  (  <{>^-(j)^)  +  <j>gcosi|jccos  (  <!>G- 4>2  ) 

-COS  (  <\>  -<f>  )  -  iji^COS  ij)  sin(  —  6-  ) 

U  Z  e>  C  (j  c  2 


(3.2.4) 


r  (  2  )  , 

i  n  ^  J  - 


cos  (  /t>G-<Pc-*2) 


-sinU^c-V 


(3.2.5) 


Equation  of  Meshing 

The  sliuiny  velocity  V  is  represented  by  the  following 


equation 


V(  21)  =  J2)  x  ( j  0,  +r[2)  ) 
~  ~  2  1  ~t 


Here  (Fig.  3.1.2) 


_  (Dvr(2) 

a)  xr  ^ 


(3.2.6) 


o  ,N1+N2. 

L  =  r 2  (— nt5 


Ni 

^2  )  —  2a^  if> 


1'ne  expression  for  co^  is  obtained  by  differentiating  Eq.  (3.1.1). 
Equation  (3.2.6)  yields 


-  _YU)  - 

\ - 

o 

CJ 

1 _ 

t 

,  Ul)  _ 

V  - 

,  K  2)  UK  . 

(  co  —  0)  )  + 

0 

t 

u 

0 

(3.2.7) 


Equations  (3.1.5),  (3.2.7),  (3.2.4)  and  (3.2.5)  yield  the  following 
equation,  of  noshing 

N1  N1 

£(iG^3i)=  cos ->c  (tt—  +  1  -  2a1*1)  -  (—  +  1)  cos(  5g-^c-^2)  =  0 

x  2 

(3.2.8) 

“1  .  2 . 
wiie  re  y  -  - —  5  ^  +  ( u  i  -  a  j.  4>  ^ )  . 

“z 

rixiiij  y^,  we  can  obtain  parameter  <*>q .  equation  (3.2.3)  provides 
two  solutions  for  £  .  Therefore  it  is  necessary  to  choose  the 
solution  corresponding  to  cne  wording  part  of  shape  1^' 

Pinion  Tooth  Shape  Equation s 

matrix  L ‘-I i f J  is  represented  oy  (Fig.  3.1.2) 


lMlf  1 


costj  -sinj^  U  0 

sin-^  cos)j  6  0 

U  0  10 

U  0  0  1 


(3.2.9) 
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Equations  (3.2.9),  (3.1.7),  (3.2.4)  and  (3.2.8)  yield  the  following 

equations  for  the  pinion  shape 


x ,  =  -C°sina  +  r„sin  [a, -il»  +g  (  *  .  )  ] 


1 


1  2  1  n  1  v  c  ^  1 


(3.2.10) 


1  2 

+  r?  [—  +  b^  -  a1*1  +  il>c  -  g(d>1)]  cos  ii^cos  [  ( <b  q+g  (  *>  q  )  ] 

2 


y1  =  C  cos*1  -  r2cos  [  ^  1-j)c+g  (  *  ^  )  ] 


^  1  2 

+  r2  [—  +  bi  ~  a  q  i  +  4»c“  -3  ( 4  ]_ )  1  cos^oinUq+gUq)] 


He  re  : 


2a 

■J  (  4  L )  =  arc  cosf[l-  (~N  y5'2'  )TT]  cos"-’c  1 

Equations  (3.2.10)  represent  the  pinion  shape  in  parametric  form,  with 
the  parameter  4q.  For  the  case  where  the  prescribed  transmission 
errors  are  zero,  a^  =  0 ,  b  q  =  0,  g  ( $  q )  =  i'i  and  Eq .  (3.2.10)  represent 
the  pinion  tooth  shape  as  a  regular  involute  shape. 


Pinion  Tooth  Surface  Equations 

The  crowned  pinion  tooth  surface  must  be  deviated  from  a 
cylindrical  surface  in  the  longitidunal  direction  (Fig.  1.3).  The 
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simplest  way  to  outain  such  deviation  is  based  on  representation  of 
the  pinion  tooth  surface  as  a  surface  of  revolution.  Such  a  surface 
can  be  _,onoruteu  oy  tne  rotation  of  the  pinion  tooth  shape  about  a 
tixeu  axis.  The  determination  of  the  pinion  surface  of  revolution  can 
oe  performeu  oy  using  tne  following  procedure: 

(i)  Assume  that  the  pinion  tootn  shape  has  been  transmitted  from 
tne  coordinate  system  Sn  to  the  auxiliary  coordinate  system  Sa  (Fig. 
3 . .  la )  .  Thus  the  tootn  shape  is  now  rigidly  connected  to  coordinate 
system  Sa  and  is  represented  in  Sa  by  the  following  matrix  equation: 


lr  j  =  [ M  .  ]  [ r  (  6 1  )  ] 
a  a  1  11 


(3.2.11) 


wne  re 


LlialJ  = 


(3.2.12) 


(  li j  on  tne  second  sta^e  we  use  a  fixed  coordi  ite  3^  ana  rotate 
tne  coorainate  system  Sa  witn  tne  tooth  shape  about  the  yD-axis.  The 
toouh  snape  wiii  generate  a  surface  of  revolution  represented  in  as 
l  jliowi)  (t  1j.  3  .  /. .  1  o  ) 


ir  j  =  [M  J  ir 

ij  O  d  ci 


(3.2.13) 


.  '.V  V. 

y*  /»///*/' S' y*. 


-  -  v- *  %\-Y A V-  .*->VvvV 
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[M  J  = 
oa 


cos  9 
0 

sine 

U 


-sme 

0 

COS  9 
0 


(3.2.14) 


(lii)  bow,  we  can  represent  the  generated  pinion  tooth  surface 
coordinate  system  using  the  following  matrix  equation  (Fig. 


3 . zic  i 


=  lHlo^£ 


(3.2.15) 


wnere 


LMlbJ  = 


(3.2.16) 


Tne  final  expression  of  the  pinion  tooth  surface  can  be 
represented  now  as  follows 


Li'i  (  9  ,  <t> .  )  j  = 


f x( 9, *  ) 

f .  (0,4>..) 

^  p 

f  3(e,y 


[:ilb]lWba(9)j[Maljtrl(^P)J 


(3.2.17) 


liquations  (3.2.  11)  -  (3.2.17)  yield 


f  i  (  9  ,  <}>  )  =  xx(  $p)cos  9 


R( 1-cos  9 ) 


=  yxU y) 


£  j  (  9  ,  <t>  p )  =  x^l^Jsine  +  Rs  i  n  9 


(3.2.1 8) 


here  x .  ( <p  ),  /,(<+)  )  are  functions  (Eg.  3.2.10)  tnat  represent  the 
i  r  i  p 

pinion  shape  in  the  middle  cross-section.  Also  for  convenience,  we 
use  h>  as  curve  parameter  instead  of  <j>^. 


3.3  Simulation  of  Meshing  and  Bearing  Contact 

The  authors  have  developed  the  TCA  (Tooth  Contact  Analysis) 
program  to  simulate  the  meshing  and  bearing  contact  of  misaligned 
gears  (see  the  Appendix).  The  results  of  the  TCA  program  confirm  that 
the  proposed  method  of  crowning  reduces  the  sensitivity  of  gears  to 
cue  misalignment  (as  shown  in  numerical  example). 

Example 


oiven:  numbers  or  teeth:  N ^  =  20,  N 2  =  40;  diametral  pitch 
f  =  lu  —  ;  pressure  angle  tii  =  20°.  The  pinion  tooth  surface  has  been 
designee  as  a  crowned  surface  and  the  function  of  transmission  errors 
lor  tne  years  without  misalignments  has  been  represented  as  a 
paracolic  function  with  the  "Level  of  Kinematical  Errors  (LKE)"  =  2 

arc  second.  Radius  of  rotation  Rj  =  350  in  (see  Fig.  3.2.1).  The 

enveloped  TCa  program  has  been  applied  for  the  evaluation  of  trans¬ 
mission  errors  for  the  following  misalignments: 

Ac 

(1)  Tno  change  of  the  center  distance  is  —  =  1%.  The  gear  axes  are 

not  paraiiel  but  crossed  and  the  twist  angle  is  5  arc  minutes. 


2u 


P’.VA  V.V  -.'AA'  V  VV  -p  -V  V  V 


-,V  V. 


The 


function  of  kinematic  errors  caused  by  the  misalignments  mentioned 
above  is  of  a  parabolic  type  and  the  maximum  value  of  transmission 
errors  is  1.2  arc  seconds. 

(ii)  The  ‘_,ear  axes  are  not  parallel  but  intersected  and  form  the 
angle  a  =  5  arc  minutes.  The  function  of  kinematic  errors  is  of  a 
parabolic  ty^e  and  tne  maximum  value  of  transmission  errors  is  2.0  arc 
seconds . 
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4. 


Deviation  of  Pinion  Tooth  Surface  by  Tool  with  Cone  or 
Revolute  Surface 


4.1  Introduction 

A  new  raetnod  for  generatation  of  tne  pinion  of  spur  gearing  is 
developed.  Tne  involute  spur  gears  with  the  crowned  pinion  tooth 
surface  are  less  sensitive  to  tne  gear  misalignment,  the  bearing 
contact  has  a  favorable  location  and  tne  kinematical  errors  are  very 
small  and  witn  favorable  snape .  The  crowning  may  be  done  by  grinding 
or  snaving  using  a  tool  with  conic  surface,  a  surface  of  revolution  or 
a  planar  surface.  The  new  method  covers:  (i)  the  theory  of  the 

proposed  method  of  crowning,  (ii)  determination  of  the  pinion  surface 
generated  oy  tne  tool,  (iii)  conditions  of  non-undercutting,  (iv) 
principal  curvatures  and  directions  of  the  pinion  tooth  surface,  (v) 
dimensions  and  orientation  of  the  instantaneous  contact  ellipse, 
(vi)  the  bearing  contact  and  (vii)  the  kinematical  errors  due  to 
.jear  misalignment.  Computer  programs  are  developed  in  all  stages 
including  the  simulation  of  meshing  and  bearing  contact  TCA  (Tooth 
Contact  Analysis)  program. 


<* .  2  Principle  of  (Generation  and  Coordinate  Systems 

consider  two  rigidly  connect  generating  surfaces  Eq  and  E  .  The 
■generating  surface  is  a  plane  and  generates  the  gear  tooth  surface 


£.  tnut  is  a  regular  involute  surface.  Surface  E  is  a  surface  of 
i  P 

revolution.  Initially  we  consider  that  E^  is  a  cone  surface  and 


l  and 


contact  each  other  along  a  straight  line  that  is  the 


jo  no  tutrix  of  tne  cone. 

Figure  4.2.0  shows  the  generating  surfaces  and  E  .  Figure  1.4 


snows  a  sKetch  of  the  tool  chat  could  oe  used  for  pinion  generation. 

Figure  4.2.1  illustrates  the  process  of  generation.  While  the  rigidly 

connected  generating  surfaces  perform  a  translational  motion,  the 

pinion  and  gear  rotate  about  their  axes  0^  and  02  respectively.  The 

parameter  of  motion  of  the  cutters,  s,  and  the  angles  of  rotation  of 

the  pinion  and  the  gear,  $  and  ,  are  related  as  follows 

p  o 


3  =  rl*p  =  r2*G 


(4.2.1) 


wnere  rj_  and  r2  are  the  gear  centrodes  radii;  the  cutter  centrode  is 
the  straight  line  that  is  tangent  to  the  gear  centrodes.  Point  I  is 
the  instantaneous  center  of  rotation.  Coordinate  systems  Sp  and  S2 
are  rigidly  connected  to  the  pinion  and  the  gear  respectively.  The 
ridiaiy  connected  tool  surfaces  (see  Fig.  4.2.0)  are  represented  in 
coordinate  system  Sc;  Sf  is  the  fixed  coordinate  system. 

The  generating  surface  (a  plane)  is  covered  with  the  set  of 

contact  lines  LG2  —  the  instantaneous  lines  of  tangency  of  surfaces 
and  (Fiy.  4.2.2a).  The  location  of  these  lines  depends  on  the 
value  of  parameter  <^.  Line  LGp  is  the  line  of  tangency  of  generating 
surfaces  z ^  and  Ep.  Points  M^,  M2  . . . ,  Mn  represent  the  set  of  points 
of  intersection  of  LGp  and  Lq2 .  The  coordinate  transformation  from  Sc 
to  S2  is  represented  by  the  equation 

ir  zj  =  ln2(J  lr^G)  ]  (4.2.2) 

n  (  °  J  -  v  (  ^  ^ )  =  lVuG'9G'  '('(j)  =  0  (4.2.3) 


23 


Here:  r'  ( U’  ,  9  )  is  the  vector  function  that  represents  plane 

is  the  matrix  that  describes  the  coordinate  transformation 
from  Sc  to  £>2*  Equation  (4.2.3)  is  called  the  equation  of  meshing  of 
tne  ^,ear  cutter  and  gear  2  [Litvin,  1987].  nv  'istheunit  normal  of 
generating  surface  and  is  the  relative  velocity 

for  tne  process  of  generation.  Contact  point  Mq  of  surfaces  E and 
L  2  may  be  represented  in  coordinate  system  Sc  as  a  function  of 
parameter  d>,,  as  follows 

Cj 


£cG>  =  EcG,'ug.9g) 


F2(UG'9C)  =  0 


f 2  ^  UG '  0G ' *G  ^  ^ 


(4.2.4) 


bquat  ions 


Ec^^C'V' 


F2(UG,0G)  0 


(4.2.5) 


represent 

line  L(jp  in 

coordinate 

sys  tern 

5C .  Equations 

(4.2. 

4) 

represent 

trie  location 

of  points  M 

1'  ri  2 '  * 

.  .  ,  Mn  on  surface 

ZG  as 

a 

f  unction 

of  parameter 

<{>r,  .  Matrix  Eq. 

(4.2.2)  and 

Eq . 

(4.2. 

4) 

represent 

tne  sec  of  contact  points 

of  gear 

tooth  surfaces 

L  .  &  i  t  Ui  l,  i 

^  1 

in  coordinate 

system  ^2* 

These 

points  are  the 

centers 

of 

tne  instantaneous  contact  ellipses 

that  form  the  bearing 

contact 

of 

■juar  2  and  pinion  L. 

Using  the  matrix  equation 


>lr  ]  [r^> 

f  c  c 


(4.2.6) 


ana  Bq.  (4.2.4)  we  .-nay  represent  in  coordinate  system  Sf  the  line  of 

action  -  tne  set  of  contact  points  of  E ^  and  E ^  in  the  fixed 

coordinate  system. 

similarly,  considering  the  generating  surface  E  ,  we  may 

hr 

represent  in  coordinate  system  Sc  the  lines  of  contact,  Lpi'  of 
surlaces  Ep  and  E  ^  as  follows 

£cP)  =  EcP)(uP,0P)  n(P).v{Pl)  =  f1(upr0p/*p)  =  0  (4.2.7) 

(  P ) 

ne re  n  is  the  unit  normal  of  generating  surface 

E  ;  V<Pl>  =  v^P^  -  is  the  relative  velocity  for  the  process  of 

r  ~  ~  ~ 

generation . 

Tne  pinion  surface  is  represented  by  the  equations 


Lci  i  = 


lrflcJLrc 


(P) 


=  [MlfJ[HfcJ[rcP)] 


f ,  (u„,  9 


p^p 


)  =  0 


(4.2.8) 


Lqua  t ions 


.(  P) 


£c  P }  (  U  P  '  9  P  }  ' 


F(up,9p)  =  0 


(4.2.9) 


it,  in  terms  of  up,  and  0D,  the  line  of  contact  of  surfaces 


and  E.,.  in  coordinate  system  S  ,  Contact  point  Mi  of  surfaces 


E,  and  E  may  be  represented  in  coordinte  system  S  as  a  function  of 
1  z  ^ 


parameter  <f>p  as  follows 


(p) 


r.c  '  Ec'’’  (ivv 


i  i  (UP  >  °P'  4>p  )  0 


(4.2.10) 


ivuB' V 


u 


2d 


Figure  4.3.2  illustrates  the  installment  of  the  tool  cone  in 


coordinate  system  Sc .  Axes  of  coordinate  system  are  parallel  to 
axes  of  system  Sa  and  axes  x^ ,  y^  lie  in  plane  zc  =  0.  The  coordinate 
transformation  from  Sa  to  3C  is  represented  by  the  following  matrix 
equation 


1Hob‘ ‘"ta1 lra] 


"  1Hca1!ral 


(4.3.4) 


IMcb> 


COS  (at—  ^  ) 

-sin(ct-ij4  ) 


sin  ( a-t|»  ) 

cos  (  a  -  ) 

c 


(4.3.5) 


a  o 


-upsina 

* 

-  (d-UpCosot ) 


(4.3.6) 


U  0 


Comoininy  Eq.  ( 4 . 3  .  t> )  and  Eq.  (4.3.6)  results  in: 


[McaJ  = 


COS  (  Ct- lpc  ) 
■Sin  ( <x~  ) 


sin  (  ct  —  ^  ) 
cos  (  ) 


-upsimj»c-dsin(  ct- 1|»  ) 

★ 

u  pCOS  6c~dcos  (  a-  i|<  ) 


(4.3.7) 


*-■  v /  v,  /  -* 

-,vV 


•  *  .•  .  .*  /  /  /  „*  .V  W~.  W_  V.  mT  f  W  tfb  r.  w-v  « 

-.%v.vr 


Here  u„  =  I  AO 
r 


b 


cos  <ii 


Using  Eqs.  (4.3.4),  (4.3.7)  and  (4.3.1),  we  represent 
the  generating  surface  Ep  in  coordinate  system  Sc  by  equations  as 


follows 


^p,<  - 


u  p  [ s inacos  0 pcos  (  a  —  ij;  )  -  cosasin  ( at-i^c  )  ]  -  upsinij*c 

k 

-u  dsinacose  sin(a-if)  )  +  cosacos  (  a-'p  , )  ]  +  upcos<K 
u^sinasine^ 


(4.3.8) 


It  is  evident  that  the  point  cf  the  cone  surface  with  the  surface 

* 


origin 


0, 


of 


parameter  up  =  up  and  6p  =  U  coincides  with  the 
coordinate  system  Sc .  The  unit  normal  of  cone  surface  is  represented 
in  coordinate  system  Sc  by 


,n(P) 

inc  J 


iL  Jin] 
ca  a 


(4.3.9) 


ncre 


l  u  J 


cos  (rt-\ ii  ) 
-3  in  (  a  -  ) 


sin  ( a-i]/c) 
cos  (  a-i  1>  ) 


0 

0 

1 


(4.3.10) 


r 


‘n‘P,l  ■ 


cosaCOSBpCOG  (a-ti;c  )  +  sinasin(a-tbc) 
-co3aco59pSin(c«-’l)c)  +  sinacos  (a -if;  ) 
cosusing 


(4.3.11) 


4 . 4  Pinion  Tooth  Surface 

The  equation  of  meshing  of  the  generating  surface  zp  and  the 
pinion  tooth  surface  is  represented  by 


n(P>.V(P1)  =  n(P).(V(P)-V(1>)  =  0 


(4.4.1) 


We  nav  also  use  the  equation  based  on  the  condition  that  the 
contact  normal  must  intersect  the  instantaneous  axis  of  rotation  I-I 
(Litvin,  1967).  Thus  we  obtain  (Fig.  4.2.1) 


(  P) 

A  -X 

c _ c 

n(P>‘ 


Y  -y<P) 
c  1  c 


Z  —  z ' 
c  c 


(4.4.1a) 


Here:  (Xx,  Y c,  Zc)  are  the  coordinates  of  a  point  that  lies  on  axis 

(  p )  (  P )  (  P ) 

I-I;  X  ,  Y  and  Z  are  the  coordinates  of  the  cone  surface;  n„Y, 
c  c  C  ^ A 

n  and  n  are  the  projections  of  the  surface  unit  normal.  Here 
(Fig.  4.2.1) 


3  =  rl*P 


Yc  =  0 


(4.4.2: 


Lquations  (4.4.1a)  (4.4.2),  (4.3.3)  and  (4.3.11)  yield 


■ -.’tvv-vvr v\«.‘ ■ *.*  ■ sv/y> 


1 1  ( up  '  9p  '  4>p  )  =  r^pi-cosoicos^sinta-^)  +  sinacosfa-^)] 


*  2  2 
•Upcos0p  +  up( cos 0 pCos  a  +  sin  a)  =  0 


(4.4.3) 


r.  ■r 
«\- 
■ r.* 

m 


■v** 

L'v*- 

vx 

nsfr 


ft 

K'.\ 


here  (  Fi^ .  4.3.2) 


*  cl 

u  ,  =  At)  = 


a  cosa  cos  ip 


(4.4.4) 


Tne  line  of  contact  LpG  of  generating  surfaces  E D  and  z( 


(Fig.  4.2.0)  is  generatrix  of  the  cone  surface  determined  with 


9P  =  0 


(4.4.5) 


Equations  (4.3.0),  (4.4.3)  and  (4.4.5)  represent  the  location  of  the 

instantaneous  contact  point  of  surfaces  E  ^  and  E2  on  surface  Ep  as 


follows 


(  P) 


=  r 


j_  $  p5^  h  ^ 


(P)  _ 


-  r^pSin^cos^ 


z(P)  =  0 
c 


(4.4.6) 


1'ne  coordinate  transformation  in  transition  from  Sc  to  S£  is 
represented  oy  tne  matrix  equation 


it,J  =  [MfJ[r(P)] 

JL  LL. 


(4.4.7) 


wnere  (  Fi  g  .  4 . 2  .  1  ) 
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f  -*  /.  .*  -*  S.  /  v  .*  .*  V 


•  -r. 


f-/-/  - 


%  * 


0  1  u 
0  0  1 
o  u  0 


(4.4.3; 


D^uations  (4.4.0)  to  (4.4.0)  yield 


X£  ~  ri^pco£i  ->c 


Y  f  =  r1_r1  ♦pSin^cos^ 


zf  =  0 


(4.4.9) 


Equations  (4.4.9)  represent  the  line  of  action  as  a  tangent  to  the 
case  cylinder  of  radius  r^  that  lies  in  plane  Zf  =  0  and  passes 

through  point  I  (Fig.  4.4.1). 

The  pinion  tootn  surface  may  be  represented  in  coordinate  system 
S  y  as  the  set  of  contact  lines  of  surfaces  £p  and  Thus  we  obtain 


“V  =  iflicllrcP)i  =  lMi£’  lafC|  [rcP)| 


(4.4.10) 


nere  ( Fi j •  4.^.1) 


i;,ifJ 


cos  $ 
-s  i  n  it 


sinip 

cos 


(4.4.11) 


h juat ions  (4.4.10),  (4.4.3)  and  (4.3.8)  yield 


3.3V-. 


WIVUWJH 


V: 


m 


■»  .< 


m 


m-. 


m 


S&:: 


r.v, 

►»'  V 


f>X< 


h'v 


y.  ■■. 


f-v* 

a  „• 


J& 


K;I 


xp  -  Up  [cos  e^sinacos  (  a-ijic+<tp)  -  cos  as  in  (  a- ^  +  <f>p  )  ] 


-Up  s i n  (  0C ~ <!>p )  +  rx  (-<*,pcos  <j,p+sin<))p; 


y  ^  =  -  Up  [cos  6p3  in  as  in  (  a- '!>c  +  <J>p  )  +  cos  acos  (  a- +  <£p  )  J 


+  Up  cos(^c-4,p)  +  r1(^psin«p  +  cos<,p) 


-  u  sintpSina 


f  1  (Up  ,  9p,  -}>p  )  =  rl  [_c°sacos  0psin  (  a-tjj  )  +  sinacos  (  a-  )  } 


*  2  2  (4412) 

-u„cos9^  +  un(cosencos  a  +  sin  a)  =  0 

r  r  r  r 


equations  (4.4.12)  are  equivalent  to  the  representation  of  the  pinion 
suriace  by  the  vector  function  r  p  (up,  9p, <j> p )  and  the  equation 
of  raesning  f  ^  ( u  ,  3p ,  <f>p )  =  0.  Taking  it  into  account  that 


l  1  ( j  (  B  p /  i> p )  -  u  is  linear  with  respect  to  the  parameter  up,  it  is 


easy  to  eliminate  up  and  represent  the  pinion  surface  in  two 
parametric  form,  with  tne  parameters,  9p  and  <f>p. 

me  intersection  of  tne  surface  by  the  plane  zp  =  0  represents  a 
regular  involute  curve.  aquations  of  this  curve  may  be  derived  from 
Lq.  (4.4.12),  taking  it  into  account  that 


i  U  b  oi  L  L  « 


9  .  =  u 

t 


ana 


Up  u  p 


r  p  it>  pS  1  n  ih  ^ 


(4.4.13) 


I  n u  s ,  wo  nave 


x,  =  r  I-.}  cos4  cos(  )  -!],  )  +  sin*  ; 

line  r  C  r 

/  !  =  pcos  ycs  i  n  (  yp-  ^  )  +  cos  xp] 


(4.4.14) 
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am »*■> wV VN.  .>  »y  a *-VkA  J~i  -V ^ -■ 


const,  we  obtain  the 


Cutting  the  pinion  tooth  surface  by  a  plane  z-^  = 
cross-section  of  the  pinion  tooth  surface  that  is  deviated  from  a 


regular 

involute  curve 

(Fig  . 

4.4.2) . 

It 

will  oe  shown 

iater 

that  the 

kinematic  errors 

caused 

by 

tne 

>j  (j  ci  r  rn  l 

sali^nment  can 

be  of 

two  types 

that  are  shown  in 

Fig .  1 

•  2a 

and 

L.2u.  inertia  tic  errors  shown  in  Fig.  1.2a  correspond  to  the  case  when 
tne  tranof  ormation  of  rotation  is  interrupted  and  the  gear  tooth 
surfaces  are  out  of  contact.  To  avoid  the  appearance  of  such 
ft i nematic  errors,  a  surface  of  revolution  that  slightly  deviates  from 
tne  cone  surface  is  actually  used.  The  determination  of  curvatures  of 
tne  surface  of  revolution  is  a  subject  of  optimization.  Since  a 
surface  of  revolution  is  used  instead  of  a  cone  surface,  the  crowned 
pinion  surface  deviates  from  a  regular  involute  surface  in  all  cross- 
sections  including  z =  0,  the  middle  cross-section  of  the  tooth 
surface  . 


4.0.  Conditions  of  Pinion  Without  Undercutting 
tencrai  Approach 

Tne  problem  of  undercutting  of  the  pinion  tooth  surface  by 
crowning  is  related  with  the  appearance  on  the  pinion  tooth  surface  of 
siiijiiiur  joints.  Lt  is  Known  from  differential  geometry  that  the 
suL'j-aeo  point,  is  singular  if  the  surface  normal  is  equal  to  zero  at 
suc.i  u  ljU  mt.  Litvin,  proposed  a  method  (Litvin,  19o7)  to  determine  a 
line  on  tne  tool  surface  which  will  yenerate  singular  points  on  the 
s  jri.au  e  that  is  jenerated  by  the  tool.  This  line  designated  by  L 
(tij.  4.1.1;  must  be  out  of  the  working  part  of  the  tool  surface  to 
•jvo  in  unuercuttinj  of  tne  pinion  by  crowning. 
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The  limiting  L  of  the  tool  surface  is  determined  by  the  following 


equations 

r ( ^ )  =  r(P) ( u  u / 0  D )  (4.5.1) 

~C  ir  r 


is  to  be  represented  in  coordinate  system 


Derivation  of  Equation 

( cl ) 

Relative  velocity 
oc  as: 


;  ( c  1 ) 


=  V(C)  -  V(1) 


(4.5.5) 


:c) 


(i)  ... 


wnore  V'  is  the  velocity  of  the  cutter  and  V  is  the  velocity  of 
~c  J  ~c 

the  pinion  at  generating  point.  Using  the  sketch  in  Fig.  4.2.1,  we 
yet 


(4.5.6) 


•  ds 
dt 

riwp 

v(c)  = 

~c 

c: 

= 

0 

i 

c 

I 

O 

1 - 

V 


(1)  _ 


=  “p  +  0  U.xwo 


(4.5.7) 


deriving  Eg.  (-1.5.7)  we  suDstituted  the  sliding  vector  up  that  passes 

through  u ^  oy  the  equal  vector  tnat  passes  through  0C  and  the  vector- 

moment  0  o,  xu,, . 

e  1  ~P 

here  : 


u 

ri*P  ' 

= 

0 

w  0  0,- 

c  1 

'rl 

1 

0 

Equations  from  (4.5.6)  to  (4.5.8)  yield 
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-  v.v.:<a 


■V 

,\5 

ti 


I 


■x '  '+r 


iyP 


(4.5.9) 


(  p  )  (  p  ) 

wnere  x  „  anu  y^  are  represented  by  equations  ( 4 . 3 .  ii ) 
(4.u.u),  (4.3.5),  (4.4.3)  and  (4.5.4)  yield 


Equations 


r  i  (  u  p ,  3 , . )  =  m^r^q  -  (m2+m  )qp  -  myn^q+m^p  =  0 


(4.5.10) 


Mere 


u  ,-,COS  Ot 
P 

___ 

* 

UpCOSij)  ^ 

d/cosa 

.a^  =  sinasm  (  a- !|;  )  cos  9p+  cosacos(a-fc) 

-  s  inaoos  (  a- -p  )  -  cos ats i n  (  a- )  cos  J  9 p 

3 

f.\  ^  -  Lsinoicos  (  a-ii»c  )  -  cosasin  (  a— tp  )  cos  3p] 

.  2 

r  a  =  s ina cosasin  3 

4  P 

r  ^  cos  at 

5  =  __J 

qFl 

using  tq .  (4.5.10)  and  taxing  it  into  account  that  ■  *  0,  we 

a  9p 


4  “ 


P  = 


d,(9J  >  — 

P  P  cosa 


(4.5.11) 


ine  analysis  of  Lq.  (a.5.1U)  yields  that  inequality  (4.5.11)  is 
satisfied  if  tne  conditions  of  non-undercutting  of  the  pinion 
generated  by  a  reyular  rack  cutter  are  satisfied. 

4.6.  Principal  Directions  and  Curvatures  of  Tooth  Surfaces 
Introduction  to  Pinion  Principal  Directions  and  Curvatures 

The  direct  determination  of  pinion  principal  directions  and 
curvatures  requires  complicated  derivations.  A  simplified  approach 
for  the  solution  to  this  problem  has  been  proposed  by  Litvin,  in 
livj.  The  main  idea  is  as  follows: 

Consider  that  the  tool  surface  Ep  generates  surface  E ^  as  the 
envelope  of  the  family  of  surfaces  Ep.  Let  Lp^  be  the  instantaneous 
line  of  cuntact  of  Ep  and  E  ^  and  M  is  a  point  of  this  line.  We 
consiuer  as  given  tne  principal  curvatures  and  directions  of  Ep  at  any 
point  and  also  the  parameters  of  motion  of  surfaces  Ep  and  E ^  being  in 
mean.  Then  tne  principal  curvatures  and  directions  of  the  generated 
surtace  E  ^  can  be  expressed  in  terms  of  principal  curvatures  and 
directions  of  tne  generating  surface  Ep  by  using  the  equations  that 
nave  seen  derived  by  Litvin,  1969. 


Principal  Curvatures  and  Directions  of  E 

- & - - - - — -  p 

The  tool  surface  is  a  cone  surface  and  its  principal  directions 
cai no i ue  witn  tne  direction  of  the  cone  generatrix  and  the  direction 
mat  is  perpendicular  to  the  cone  generatrix.  The  Rodrigues'  formula 
(  r:.g .  4.0.1)  can  do  used  to  find  the  principal  curvatures  and 


(4.6.1) 


directions . 


Here:  <T  rT  are  the  principal  curvatures  of  surface.  V  is  the  velo- 
i,ii  ~r 

city  of  a  point  that  moves  over  a  surface  and  rt  is  the  derivative  of 
the  surface  unit  normal  n,  when  n  changes  its  direction  due  to  the 
motion  over  the  surface. 

equation  (4.u.l)  yields  the  following  expressions  for  the 
principal  curvatures  and  directions  (see  Eq.  (4.3.1)  for  the  tool  cone 
surface  represented  in  coordinate  system  3  ) 

tr  a 


(P) 

i 


-sin  9p 
0 

cosep 


K^)  = _ I_ 

I  uptana 


e 


(P) 

II 


cos0pSin« 

-cosa 

sinepsina 


(4.6.2) 


(4.6.3) 


(4.6.4) 


(4.6.5) 


i'ne  negative  sidn  for  indicates  that  the  curvature  center  is 

located  on  tne  negative  direction  of  the  surface  normal. 

fne  principal  curvatures  are  invariants  with  respect  to  the  used 


coordinate  system.  The  unit  vectors  of  the  principal  directions, 

(  P  )  (  P ) 

e'  arad  e'  ,  may  be  represented  in  coordinate  system  Sf  by  using  the 
following  matrix  equation 


i  t  X.  i  Ju 


=  [Lc  ] [ L  ] [e$P|_ 
fc  ca  I ,  Ii 


(4.6.6) 


(  P )  (  P ) 

here:  ie'  'j  and  [e r  iT]P  are  vectors  represented  in  S_  and  Sf, 

i ,  x  x  ci  1,11  r 

respectively;  matrix  [Lca]  is  represented  by  Eq.  (4.3.10);  matrix 
Lt£cJ  is  tne  3x3  unitary  matrix  (see  matrix  Eq.  (4.4.8)). 

Equation  (4.6.6)  yields 


-sine  cos  (  a~i}>  ) 

ir  O 

sinopsin  ( oc  —  ) 

COS0n 


(4.6.7) 


[e(P)l 
e  1 1  Jf 


cosQpSinacos  (  a-^c)  -  cosasin (  a- \pc  ) 
-cos9psinasin(a-ij>c)  -  cosacos  ( a-#  ) 


sinGpSina 


(4.6.8) 


Principal  Curvatures  and  Directions  of  Pinion  Tooth  Surface 

Tne  determination  of  principal  curvatures  and  directions  for  the 
pinion  tootn  surface  is  based  on  the  following  equations  (see  Litvin, 
1  v  o  7  ) 
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2d13D23 


(4.6.9; 
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3cos  2cr 
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(4.6.10) 


( 1)  ,  (1) 

KII  +  KI 


(P)  .  (P)  ,  b 1 3+b  23 

1  11  b33 


(4.5.11) 


(  P  )  (  p  ) 

Flore:  and  arc  the  principal  curvatures  of  the  tool  surface 

(  P  )  (  P ) 

En;  e'  and  o'  are  the  unit  vectors  of  the  pnncioal  directions  or- 

r  ~  1  'vll 


Ip  (Fig.  4.6.1);  and  *11^'  and  eii^  are  the  principal 

( Pi  ) 

curvatures  and  directions  on  the  pinion  tooth  surface  and  a  is 

[D]  (1  ) 

measured  counterclockwise  from  e'  to  e '  (Fig.  4.6.1).  The 
coefficients  b^,  b23  and  bgg  have  been  derived  from  equations 

represented  in  [Litvin  1937]  but  modified  for  the  case  when  a  rack 
cutter  generates  a  gear.  The  expressions  for  b-^,  b23  and  b-3-3  are  as 
f o 1  lows 


c  (  P ) 

Si i *  ~p 


_0(p) 

Si  *“o 


(4.6.12) 
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(4.6.13) 


1  1  P 


ctors  that  are  used  in  Los.  (4.6.12)  and  (4.6.13)  are 


represented  in 


coordinate  system  Sf  (Fig.  4.2.1);  i  ^ ,  j  ^  and  k^  are 


tne  unit  vectors  of  coordinate  axes  of  system  Sf .  The  expressions  of 
tnese  vectors  are  as  follows 


u)  -  ok 


(4.6.14) 


is  tne  angular  velocity  of  the  pinion  beiny  in  mesh  with  the  rack 
cutter 


V(P)  = 
tr 


"urlif 


(4.6.15) 


is  tne  transfer  velocity  of  a  point  of  the  rack  cutter  that  performs 
translational  motion;  r^  is  the  radius  of  the  pinion  centrode 
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the  pinion, 

►  The 

term , 

"transfer  veloc 
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tooth  surface  is  considered 
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point 
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the  motion 

with 

the  tootn  surface 
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iiiurnj"  velocity — tne  velocity  of  a  point  of  the  rack  cutter 
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with  respect  to  tue  same  point  of  the  pinion.  Vector  V  can  be 
represented  in  terms  of  components  xc  and  yc  using  matrix  Eq.  (4.4.8) 
tnat  describes  tne  coordinate  transformation  from  Sc  to  Sf.  Then  we 
obtain 


xf  =  xc  "  rl*P 


y  +  r , 

2  c  1 


(4.6.13) 


Equations  (4.0.17)  and  (4.6.18)  yield 


V(P1)  =  „  -xc+rllp 


(4.6.19) 


Tne  principal  curvatures  oC  the  cone  that  is  rigidly  connected  to  the 
racx.  cutter  are  represented  as 


u  .,tanci 
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(P)  _  ,, 

<H  -  J 


(4.6.20) 


l’ne  unit 


vectors  of  principal  directions  are  designated  by 


(  p )  ( p )  ( PI )  ( Pi ) 

e^  an..  ej.  .  Components  Vj.  and  Vpp  are  represented  as  follows 


V.(P1)  =  v(Pl).e(P) 
i  ~  ~  I 


V<Pl>  =  V(Pl).e<P> 
11  ~  ~  1  J- 


(4.6.21) 


Alter  substitution  we  ootain 


jsin9p  up  2  rx6p 

o  ,  =  — [1  -  —  cos  a  +  -  cosas  in  (  a- '!;  )] 

1j>  s  i n ct  a  p  Up  c 


(4.6.22) 
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l  =  uPsin  9pCosasinoi  -  upsinacosa  (  1-cos  0p ) 


(4.6.25) 


k  2  =  r  i  6p  lcos  9pc°sacos  (  a  —  ^  )  +  s  i  nets  i  n  ( ct- i)j  )] 


(4.6.26) 


k3  =  r  3  icos0pcosasin(  a-i>c)  ~  sinacos  (  ct-i|>  )] 


(4.6.27) 
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(4.6.28) 


tor  tne  contact  point  located  in  tne  middle  section  of  the  pinion  we 
nave  9D  =  6  and 

u  b23  =  a)  b  3  3  =  w2  [r^pcos.j^-r  j.siniji  ]  (4.6.29) 

inen  we  obtain  the  following  expressions  for  the  pinion  principal 
curvatures  and  directions 
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(4.6.30) 


ine  principle  curvature  k^'  is  the  same  as  the  principal  curvature 
ot  a  regular  involute  surface.  However  ^  differs  from  zero 

oecause  the  pinion  is  generated  by  a  cone  but  not  a  plane  as  in  the 
case  of  a  regular  involute  surface. 


Gear  Principal  Curvatures  and  Directions 

Tne  determination  of  the  gear  principal  curvatures  is  based  on 
tno  Leiation  oetween  the  principal  curvatures  for  spur  involute 
j  e  a  r  s .  Figure  s.o.z  shows  that  two  mating  involute  shapes  are  in 
taiijcacg  at  guint  d.  aet  us  designate  the  principal  curvature  for  the 
pinion  og  <  ana  < | | ^  where  (see  Eq.  4.6.29) 


J  1)  =  _ _ K(  1)  = _ 1 _ 

I  utani>c  II  r  ^sini^-r  ^  <f>pCos  iji 


(4.6.31) 


ine  negative  sign  of  and  means  that  the  radii  of  curvatures 

are  opposite  to  the  direction  of  the  surface  normal.  The  principal 

joo r  curvature  =  0.  The  radii  of  curvature  of  spur  involute  gears 

are  related  as  follows  (Fig.  4.6.2)  where  =  -  - . - - — 

d  II  r.sinu  -r,  *ncos  i'i 

1  c  IP  c 
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Then  we  obtain 
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(4.6.33) 


The  angle 


(21) 


between  the  principal  curvatures  is  equal  to  zero. 


4.7.  Contact  Ellipse  and  Bearing  Contact 

The  dimensions  and  orientation  of  the  instantaneous  contact 

ellipse  can  be  determined  on  the  basis  of  equations  that  have  been 

developed  by  Litvin  (Litvin  et  al.  1982,  and  Litvin  1987).  The  input 

data  for  the  computation  is:  (i=l,2)  and  e  where 

and  <  j  ^  are  the  principal  curvatures,  a  ^ 2  ^  is  the  angle  formed 

by  the  unit  vectors  of  principal  directions,  and  and  e  is 

the  elastic  deformation  of  the  contacting  surfaces  at  contact  point. 

The  bearing  contact  is  formed  by  the  set  of  instantaneous  contact 

ellipses  that  move  over  the  gear  tooth  surfaces  in  the  process  of 

meshing.  In  our  case,  when  the  crowned  pinion  is  meshing  with  a  gear 

without  misalignment,  c|^  =  e  j.  2  ^  and  =  0  (Fig.  4.7.1).  The 

axes  of  the  contact  ellipse  are  directed  along  the  n~axis  and 

•-axis  ,  respectively.  Generally  the  orientation  of  the  contact 

ellipse  i  3  determined  with  the  angle  .  In  the  case  where 
(  1  > ) 

n  -9,  is  equal  to  zero.  The  equations  to  be  used  for  the 

c  ;  ~  p  ,  tali  mi  are  as  follows 


1,(1)  (2) 

—  [  v  —  »<• 


(4.7.1) 
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(4.7.3) 


_£  |  1/2 
I  A  I 


_e  1/2 
B 
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Tne  equations  for  rj.1)  and  (i=l*2,)  have  been  represented  in 

Section  4.U.  Figure  4.7.2  illustrates  the  bearing  contact  for  spur 
jOars  with  tne  crowned  pinion  tooth  surface. 


4 .  u .  Simulation  of  Meshing  and  Determination  of  Kinematical  Errors 
Coordinate  Systems 

Tne  simulation  ol  meaning  is  a  part  of  tne  computer  aided  TCA 
program  developed  oy  tne  authors .  The  simulation  of  meshing  is  based 
on  equations  tnat  provide  the  continuous  tangency  of  contacting 

O  -A  r  L  cl  J  ( i  s  • 

uenceiortn  we  will  use  tne  following  coordinate  systems:  (i)  Sf 

tnat  is  rigidly  connected  to  the  frame  and  (ii)  an  auxiliary 

c dominate  system  Sn  tnat  is  also  rigidly  connected  to  the  frame,  and 

(ill)  systems  dj  and  02  that  are  rigidly  connected  to  the  pinion  and 

tne  go  a  r ,  respectively.  The  origins  and  coordinates  axes  and  of 
coordinate  systems  dj  and  S  ^  coincide  with  each  other  (Fig.  4.3.1a). 


'-ry.y.y- vy.p-.poy.- v 


Axis  Zf  is  the  axis  of  rotation  of  the  pinion 


Also,  the  origins  and 


coordinate  axes  z^  and  Z2  coincide  with  each  other  and  is  the  axis 
of  gear  rotation  (Fig.  4.3.10).  The  errors  of  assembly  of  the  gears 
are  simulated  with  the  orientation  and  location  of  coordinate  system 
bfl  witn  respect  to  S^.  Figure  4.3.2  shows  the  orientation 
(a)  Tne  gear  axes  are  crossed  (Fig.  4.3.2a) 
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(4.3.1) 


(o)  The  gear  axes  are  intersected  (Fig.  4.3.2b) 
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(4.8.2) 


Pinion  Tootn  Surface  Equations 

Ine  pinion  tootn  surface  has  been  represented  in  coordinate 
system  by  Fq.  (4. 4.  12). The  coordinate  transformation  form  S-^  to  Sg 
is  represented  as  follows 
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LM 


f  1 


[  r  1 J 


(4.8.3) 
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Si 
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(4.0.4: 


nero  ^  is  tne  angle  of  rotation  of  the  pinion  being  in  mesh  with  the 


q  e  a  r  . 


Equation  (4.4.11)  represent  the  pinion  tooth  surface  in  a  three- 
parametric  form  in  terms  of  parameters  up,  9p  and  (j>p  where  tj>p  is  the 
angle  of  rotation  of  the  pinion  beiny  in  mesn  with  the  coned  cutter. 
Tnose  parameters  are  related  by  the  equation  of  meshing  and  only  two 
of  them  are  independent.  It  is  easy  to  eliminate  parameter  up  and 
represent  the  pinion  tootn  surface  in  S-^  in  the  two-parametric  form. 


~  1  £i('(,p'9p] 


(4.8.5) 


wnere  $p  and  9p  serve  as  the  surface  parameters.  The  pinion  tooth 
surface  unit  normal  can  oe  represented  in  coordinate  system  as 


=  j_J  [l^  J 


(4.8.6) 


.no  jxj  i.uu'ix  [  u  ij 


^  I  i  can  be  developed  from  oy  eliminating  the 


Last  row  and  tne  Last  column.  The  surface  unit  normal,  n^,  is 
represented  in  oy  tne  following  matrix  equation 


i  n  j  n 
I  c  c 


'Llt"4cl("cP)i 
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Here:  rv  is  represented  by  Eqs .  (4.3.11),  [Lfc]  is  the  3x3  unitary 

matrix  and  [L-^]  is  the  3x3  submatrix  of  matrix  given  by  Eg. 

(4.4.11). 

Tne  final  expression  for  n^^  is  as  follows 

[  n  ^ 1 } J  =  [Lfll[Llf][Lfc][n‘P)]  (4.3.8) 

Here:  [L£cj  ■'■s  tne  3x3  unitary  matrix  [I],  but  [L^]  [L^f]  t  (I] 

because  the  elements  of  these  matrices  are  expressed  in  different 
terms  ~  and  <p  ^ 

Gear  Tooth  Surface  Equations 

Tne  year  tooth  surface  is  a  regular  involute  tooth  surface  and 
may  oe  represented  in  coordinate  system  S2  by  equations  that  are 
similar  to  Eq.  (4.4.14) 

r2  i-^cosi^cos  (  <j>G"i>c  )  +  sin^G] 

r2  [ 'i>Gcos  iicsin  (  ^  +  cos*g^ 

z?  =  ).  (4.3.9) 

wnere  <>G  is  tne  angle  of  rotation  of  the  gear  being  in  mesh  with  the 
racK  cutter.  Tne  gear  tooth  surface  unit  normal  is  represented  by  the 
equations 
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3  r  ,  3  r 

x  3 ~ 


(4.3.10) 


equations  (4.3.9)  a  no  (<i.3.1u)  yield 


I' 


(4.3.11) 


To  represent  tne  gear  tooth  surface  and  its  unit  normal  in 


ItiH 


coordinate  system  Sh  we  use  the  following  matrix  equations  (Fig. 


4.0.2) 


<XB d 


^n2)j  =  [Hh2J  [r2] 


L  d  ,,  J  = 
n  z. 


s  1  n  <j) 


'"h2’1  -  lLh21[n2' 


(4.8.12) 


costj>2  -sin$.^  U  0 


COSiJ. 


(4.8.13) 


rJ> 


m 


wnere  $  ,  is  tne  anyie  of  rotation  of  the  gear  being  in  mesn  with  the 
pinion.  Equations  (4.o.b),  (4.8.11),  (4.8.12)  and  (4.3.13)  yield 


r2  {4>Gcosij>  cos  i  (  <J>2-<J>G)  +  '!<c  )  J  +  sin  (  *2- <j>G)  } 


r2  { 4>Gcos  ^sin  [  (  4>2--}>G  )  +  )  ]  -  cos  (  <b2~  *G)  } 


7J  =  \ 

LI 


1  nn  J 


-cos  [  (  <j>2~  (j>G)  + 

-sin  [  (  f  )  +  b  J 

z  u  c 


(4.3.14) 


( (4.3.  Id) 


Those  equations  with  parameters  §2  ~  6C  =  Of  -  X  =  0  determine 

point  H  with  coordinates  (0,  r 2>  0)  and  the  surface  unit  normal  at  M 

with  components  (-cos^,  -  sin^c,  0)  (Fig.  4.8.3). 

To  determine  the  year  tooth  surface  equations  and  the  surface 
unit  normal  in  coordinate  system  Sf  we  use  the  following  matrix 
equations 


(  2) 
Lrf 


“  lHfiv 


[nf2)j  =  [Lfh][nh2)j 


(4.8.16) 


Matrix  iM^j  (see  Eqs .  (4.8.1)  and  (4.8.2))  describes  the  coordinate 
transformation  from  to  for  the  case  where  the  gear  is  misaligned 
witn  respect  to  the  pinion.  Figure  %.6.4  illustrates  the  case  where 
trie  years  are  not  misaligned,  the  shortest  distance  of  gear  axes 
C  =  rx  +  r2  and  the  tooth  surfaces  of  the  pinion  and  the  gear  contact 
each  other  at  point  M.  The  misalignment  of  the  gear  will  not  cause 
tne  surface  contact  at  the  tooth  edge  since  the  pinion  tooth  surface 
is  crowned  and  deviates  from  a  regular  involute  surface. 


Equation  of  Tanqency 

The  contact  of  gear  tooth  surfaces  is  simulated  in  the  developed 
TCa  program  by  the  following  equations 


£q  1 }  (  <j>p/  Qpf  6^  ) 


(  <t>.,,  £,$.,) 


r4 l>  (  V  V  V 


Uf2)(*G^2) 


(4.8.17) 


(4.8.18; 
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contact  ^oint  vector  Eq .  (4.6.17)  provides  tne  equality  of 


.■  v  v  V  -■ 
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position  vectors  and  vector  Eq.  (4.8.18)  provides  the  equality  of 
surface  unit  normals.  Vector  Eq.  (4.8.17)  provides  three  independent 
scalar  equations  that  relate  4,  p ,  9p  ,  <J>  ^  ,  <|>q,  X  and  <j>2*  However,  vector 
Eq.  (4.8.18)  provides  only  two  independent  scalar  equations  since 
n  ^  ^  ^  and  n^2^  are  unit  vectors  and  |  n  ^  ^  ^  |  =  |n^^|  =  1.  Thus  Eqs. 
(4.8.17)  and  (4.8.18)  yield  five  independent  scalar  equations  as 
follows 

—  ^  ( i  =  l ,  2 ,  .  .  .  ,  5)  (4.8.19) 

Equation  system  (4.8.19)  is  the  expression  in  implicit  form  of 
functions  of  one  variable,  for  instance  of  <6^.  The  Theorem  of  Implicit 
Function  System  Existence  states: 

(i)  Consider  that  a  set  of  parameters 

P  ~  (  1  p  ,  9  p  ,  ,  <{)q  »  X  ,  *2  )  (4.8.20) 


satisfies  Eq .  (4.8.19)  and  that  the  Jacobian 

D(f  l,f  2'f  3'f  4,f  5} 

b(  ffep, 9  p ,  if q , a , if 2  ) 


(4.8.21) 


(the  partial  derivatives  are  taken  at  point  P°). 

(ii)  Then,  equation  system  (4.8.19)  can  be  solved  in  the 
neighborhood  of  point  P°  by  functions 
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X  (  d>  -1  )  eC  (4.8.22) 


fne  designation  C  means  that  the  functions  have  continuous 


derivatives  at  least  of  the  first  order.  Knowing  function  <J>2 (  ^  1  we 


can  determine  the  function  of  kinematical  errors 


A<^(  <f>x)  =  4.  ,(<t>1  )  - 


2V  T 


(4.3.23) 


4.9  Modification  of  Generating  Surface  £, 


basic  Relations 


The  proposed  method  of  crowning  provides  conjugated  year  tooth 


surfaces  and  a  localized  bearing  contact.  Using  the  TCA  program  it 


can  be  proven  that  the  kinematical  errors  caused  by  the  gear 


misalignment  are  on  a  very  low  level.  However,  the  shape  of  the 


function  of  kinematic  errors  is  unfavorable,  similar  to  the  shape  that 


nas  been  represented  in  Fig.  1.2a.  This  disadvantage  can  be  avoided 


0/  cne  modification  of  generating  surface  Xp.  The  authors  proposed  to 


use  a  surface  of  revolution  instead  of  a  cone  surface.  Figure  4.9.1 


snows  tne  axial  section  of  generating  surface  Zp  that  is  an  arc  of  the 


:ircle  of  radius  p.  Controlling  the  magnitude  of  p,  we  can  provide  the 


iavorable  snape  of  the  function  of  kinematical  errors  (Fig.  1.2b)  and 


also  control  tne  level  of  kinematical  errors.  The  solution  to  this 


problem  is  baseu  on  tne  relation  between  the  principal  curvatures  and 


directions  tor  tne  contacting  surfaces  and  tne  derivative 


dT~  ^ la 2 1  ^  1  ^  '  Tll^s  relation  (Litvin,  1937)  is  represented  as 


follows 


ryy. 


bll 

b  1 2 

b  1 3 

(4.9.1) 

b21 

b22 

b23 

=  0 

b31 

b  32 

b33 

Consider  that  the  tooth  surfaces  of  a  crowned  pinion  and  a  gear  are  in 
contact  at  point  M  that  coincides  with  the  instantaneous  center  of 
rotation  (Fig.  4.0.4).  It  is  assumed  that  the  gears  are  not 
misaligned.  In  this  case  the  coefficients  of  determinant  (4.9.1)  are 
expresseu  as  follows 
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b33  =  -Q,  (1)xv|J))  -  («o(2)xV+^1)]  +  (M(1)]2ra£1(nxk2).(r^1)-Rf 

(4.9.2) 

where 

at  -  —  -  m'  =  — -1  =  _iL  (— -) 

21  uiji^  ’  21  d  ^  d  <|>  ^  d  (j>  1 


equations  (4.3.11),  (4.0.U),  (4.0.9)  with  ep  =  0  yield  (see  Fig.  3.4) 
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equations  (4.9.1)  and  (4.9.4)  yield  tne  following  relation  between  the 
derivative  and  the  principal  curvatures 


in 


21 


^ 1  >  ^ 1  (  i  \  i  0) 

( 1+Try>  ‘  (1+hj)  sin*criI'KIi  +  nr  1 


,  U2,  ( 1)  (2) , 

rjCOS^l-Kn  +  ku  i 


(4.9.5) 
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wnero  o  is  tne  radius  of  the  arc  circle  that  generates  the  surface  of 
re 7o 1 u L ion  . 


JO 


N  S 


Tnere  are  two  important  particular  cases  for  the  application  of 


tq.  (4.S.7)  (i)  The  generating  surface  £p  is  a  cone.  Considering 

tnat  p  -»■  00  ana  the  generating  surface  becomes  a  cone,  we  obtain  that 
‘azi  =  ^  1L  tue  years  are  not  misaligned.  In  this  case,  when 

misalignment  exists,  kinematic  errors  will  occur  as  shown  in  Fig.  1.2a 
wnien  are  very  undersirable .  (ii)  Tne  generating  surface  Ep  is  a 
revolute  surface,  p  is  positive,  <  u  (  I  I  =  — )  and 

11  'll1  P 

m.'  <  u  (m'  is  taken  at  point  M)  .  Then  the  function  of  kinematic 

Z  i  Z  1 

errors  will  be  similar  to  that  as  shown  in  Fig.  1.2b  even  when  the 

gears  are  not  aligned.  It  will  be  shown  below  that  by  choosing  an 
appropriate  value  for  o  it  is  possible  to  keep  the  kinematical  errors 
ol  misaligned  years  on  a  very  low  level. 

Surface  of  Revolution  Equations 

i’ne  surface  of  revolution  is  generated  by  an  arc  of  circle  of 

rauius  p.  Ine  arc  parameters  provide  a  common  normal  for  the  cone 
sun.  ace  ana  tne  surface  of  revolution  at  tneir  point  of  tangency  M 
irij.  4.v.Zaj.  Tne  circular  arc  is  expressed  in  the  auxiliary 
cooza  mate  system  as  follows 

★ 

x^  =  p l cos ( 3 ) -cosa J  +  upsim 

*  Q,  Q, 

-  u  p  s  l  n  ot  -  2  o  s  l  n  (  a  -Hy ) 

i  -  o  i  s  i  n  (  n  +  3 )  - 

.> 

cos  ■ 


sin-xj  +  — - cosa 

Cos  <Jj 

a  c  ^ 

+  2pSilW  COS('i+-r) 
c  2 


X , ,  =  u  (4.9.9) 


Tne  installment  of  tne  generating  surface  in  coordinate  system  Sc 
is  tne  same  as  it  was  described  in  Section  4.3.  The  generating 
surface  ana  its  unit  normal  are  represented  in  S  as  follows 


l r  J  =  [M  J  [r  ]  = 
c  ca  a 


(4.9.15) 
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+  UpSinacos  (  a-<i>  )  (  cos  0p-l  ) 
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idpSina  -  2psin  —  sin(a+  -T)]sin9p 


(  p)  ,  . 

in.,  .  =  l  u  j  [n  j  = 

l,  U  a  a 


cos  9pCOS  (  a+ 3  )cos  (  a-i>c  )  +  s i n  (  a  +  B  )  s  i n  (  a- ipc  ) 
-cos  3  pCos  (  a+ 3  )  s  i  n  (  )  +  s i n  (  a  +  P  )  cos  (  a- i!i  ) 


sin  9pCos ( a  +  3 ) 


Pinion  Surface 


(4.9.16) 


Ine  ,/iiiiu:i  surface  is  represented  in  coordinate  system  as 


L  h  1  J  .-J  : 


U-,i  =  lH1£  j  iMfci  ircJ  f  (  Sj_  , 0p,  (j)p)  =  U  (4.9.17) 

sere  f(S,0p,;j>p)  -  0  is  the  equation  of  meshing  (see  Section  4.4) 
represented  as  follows 

★ 

i(S,0p<£>p)  =  r1sinficos9p  -  u  ps  i  nas  in  (  a  +  3  )  (  cos  9  -  1) 

-  r (j>psin  (  a  +  S  )  cos  (  a-’^c  )  -  cos  6 pcos  (  a  +  (1 )  s  in  (  a- ib  )  =  0 

(4.9.18) 

TCA  Program 

following  the  procedure  of  derivations  described  in  Section  4. a, 
we  can  now  uevelop  the  TCA  program  for  the  case  where  the  generating 
s a r j. a c e  Ip  is  a  surface  of  revolution.  Using  this  program,  we  can 
uei_or.il  me  tne  Kinomal  ical  errors  caused  by  the  year. 
hxdiiip  le 

b,.,  =  la  h  I  =  20,  i4  2  =  4b  C  =  l.Oltr^t^)  Ay  =  5' 

20°  a  =  JU°  o  =  500  in  d  =  0.176  in 

L're-  results  of  computation  of  kinematical  errors  are  shown  in 
r  i  j.  -+.9.  i.  It  was  found  tnat  the  function  of  kinematical  errors  has 
t  so  .  e  s  l  r  o  s  rn.ape ,  and  tne  maximal  value  is  0.16  arc  sec. 
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Conclusion 

Tne  authors  have  developed:  (i)  The  basic  principles  for  the 
optical  geometry  of  crowned  spur  years;  ( i i )  New  methods  for 
yoneration  of  crowned  pinion  tooth  surface;  (iii)  Tooth  contact 
analysis  programs  for  the  simulation  of  meshing  and  beariny  contact. 

The  features  of  the  optimal  geometry  of  the  crowned  spur  gears 
are  as  follows:  (i)  the  gear  is  provided  with  a  regular  involute 
surface;  (ii)  the  pinion  tooth  surface  is  deviated  from  a  regular 
involute  surface  to  provide  a  localized  bearing  contact  and  reduce  the 
sensitivity  of  the  years  to  their  misalignment;  (iii)  the  function  of 
transmission  errors  of  the  misaligned  years  is  of  a  parabolic  type 
witn  the  small  level  of  maximal  errors,  less  than  b  arc  seconds. 

Tne  now  methods  for  generation  of  the  crowned  pinion  surface  need 
si.upie  toois  tnat  nave  to  oe  provided  with  a  generating  plane  or  a 
SuLiaco  ol  revolution  tnat  sligntly  deviates  from  a  cone  surface.  The 
ioveiopeu  1CA  programs  provide  tne  information  on  the  bearing  contact 
and  tne  transmission  errors  for  the  misaligned  years. 
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FLOWCHART  FOR  PROGRAM  I 


START  ^ 


_ <  ’ _ 

/  INPUT  GIVEN  DATA 


r>  r>  o  o 
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*  PROGRAM  I 

*  KINEMATIC  ERROR  OF  A  REGULAR  INVOLUTE  GEAR  MESHING  * 

*  WITH  A  PINION  CROWNED  BASED  ON  PREDESIGNED 

*  KINEMATIC  ERROR  FUNCTION 

*  Vc 

*  AUTHORS:  FAYDOR  LITVIN  * 

*  JIAO  ZHANG 

*  Vc 
Vr  Vr  Vr  Vf  Vf  Vc  Vr  Vr  *  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vr  V:  Vr  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  *  Vc  Vc  Vc  Vr  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc 


PURPOSE 

THIS  PROGRAM  IS  USED  TO  CALCULATE  THE  KINEMATIC  ERROR  OF  A  CROWNED 
PINION  MESHING  WITH  A  REGULAR  INVOLUTE  GEAR.  THE  PINION  IS 
CROWNED  BASED  ON  PREDESIGNED  KINEMATIC  ERROR  FUNCTION. 

THIS  PROGRAM  IS  WRITTEN  IN  FORTRAN77 .  IT  CAN  BEE  COMPILED  BY  V 
COMPILER  IN  IBM  MAINFRAME  OR  FORTRAN  COMPILER  IN  VAX  SYSTEM. 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  2(99) , ERR (99) .ERROR (99) 

COMMON  /BLOCAL /  X ( 1 1 )  , Y (10)  , A (10 , 10)  , Y 1 (10) , IPVT ( 10) , WORK (10) , 
EPSI, DELTA, NC.NE.NDIM 

COMMON  / B L0CA2 /  ALPHA , RKS I , RP , RG , RMPG , R L , S A , CA , DR , CK , SK , HDC . C , 

S ( 3 , 3) , SF , CF , SAK , CAK . Al , XP , YP , ZP , XG , YG , R , CONST (9) 

DEFINE  PARAMETERS  USED  BY  PROGRAMS 

(1)  IN  ANG  LP  ARE  UNIT  NUMBERS  ASSIGNED  TO  THE  INPUT  AND  OUTPUT 

DEVICES 

IN-5 
LP  =  6 

(2)  NDBUG  IS  USE  TO  CONTROLL  THE  AUXILIARY  OUTPUT  FOR  DEBUGGING 
NDBUG-2 

(3)  NC  IS  THE  UPPER  LIMITATION  OF  REPEATATION  FOR  SOLVING  NONLINEAR 

EQUTIONS: 

EPSI  IS  THE  CLEARANCE  OF  FUNCTION  VALUES  WHEN  THE  FUNCTIONS 
IS  CONSIDERED  AS  SOLVED  (ALL  FUNTIONS  HAVE  FORMS  OF  F(X)-O); 
DELTA  IS  THE  RELATIVE  DIFFERENCE  FOR  TAKING  DERIVATIVES 
NC,  EPSI  AND  DELTA  MAY  BE  CHANGED  WHEN  SOLUTIONS  ARE  DIVERGENT 
OR  LESS  ACCURATE 
NC-100 
DELTA-l.D-3 
EPS  I  =  1 . D- 1 2 

(4  OTHER  P.ARAMETERS  (DON'T  CHANGE) 

NDIM-10 
NE=  5 

DR  =  DATAN  (  1  .  DO  )  /'45  .  DO 

DEFINE  INPUT  PARAMTERS  OF  PROBLEM (USE  INCH  AS  UNIT  OF  LENGTH) 


1.  PINION  AND  GEAR:  PN=DI AMETRAL  PITCH;  Nl-PINION  TOOTH  NUMBER: 
RMPG-TOOTH  NUMBER  RATIO  (GEAR  TOOTH  N0./N1); 

RKSI-PRESSURE  ANGLE;  HDC=HE IGHT  OF  DEDENDUM  OF  PINION: 
COE=CQEFF .  OF  CENTRAL  DISTANCE (USUALLY  COE-1.) 


C  R=RADIUS  OF  ROTATION  AXIS  FOR  GENERATING  PINION  SURFACE 

PN-IO.DO 
Nl=20 
RMPG-2.D0 
RKSI=20.D0*DR 
HDC- 1 . DO/PN 
COE-1. 010D0 
R-3.5D2 

C  (2)  PRE-DESIGNED  FUNCTION:  INPUT  CONSTANTS  HERE  AND  INPUT  FUNCTION 
C  IN  SUBROUTINE  FUNC 

C  MATRIX  CONST  (J)  IS  USED  FOR  CONSTANTS  OF  FUNCTION 

CONST (1) =2. D0*3. 92975 168  ID-4 

C  (3)  MISALIGNMENT:  NMIS-ID  NO. (1-CROSSING  AXES.  2= INTERSECTING  AXES) ; 
C  NG-NO.  OF  MISALIGNED  ANGLES  TO  BE  SIMULATED  (FROM  0.  TO 

C  (NG-1)*GAMMAI) ;  GAMMA I- INCREMENT  OF  MISALIGNED  ANGLE (MINUTE) ; 

NMIS-2 
NG-2 

GAMMA  I  =  5  .  DO 

C  (4)  OUTPUT:  FEEI- INCREMENT  OF  ROTATION  ANGLE  OF  P IN  ION (DEGREEE) 
FEEI-1 . 0D0,:DR 
C 

C  DESCRIPTION  OF  OUTPUT  PARAMETERS 
C 

C  FEE1-R0TATI0N  ANGLE  OF  PINION 

C  FEE2-ROTATION  ANGLE  OF  GEAR 

C  RP-RADIUS  OF  PINION  CONTACT  POINT 

C  RG-RADIUS  OF  GEAR  CONTACT  POINT 

C 

C  FIND  AUXILIARY  VALUES  FOR  CALCULATION 
RP- FLOAT (Nl)/2./PN 
RG-RP’-RMPG 
C= (RP+RG) *COE 
CK-DCOS (RKS I ) 

SK-DSIN (RKSI) 

NCOEF-360. D0*DR/FEEI/ FLOAT (N1)+0.3D0 

N=NC0EF*2+ 1 

NT-N-NCOEF 

FI  1  =  360. DO/RMPG/FLOAT (N1 ) 

C  DEFINE  MATRIX  DESCIBEING  MISALIGNMENT  AND  OUTPUT  ASSEMBLING  CONDITION 
DO  5  1=1,3 
DO  5  J=1 ,3 
S ( I , J) =0 . DO 

IF  il.EQ.J)  S ( I . J ) = 1 . DO 
5  CONTINUE 

DO  15  LL-l.NG 

GAMMA=GAMMAl*FLOAT ( LL~ 1 ) /60.D0*DR 
CG-DCOS (GAMMA) 

SG-DS IN (GAMMA) 

GAMMA -GAMMA / DR  * 6 0 . DO 
IF  (NMIS.EQ.l)  THEN 
WRITE  ILP.10)  COE. GAMMA 

10  FORMAT  (1H1.///.1X, ' C= ' ,F4.2, ' :r  (RP~RG) ;  CROSS  INC  ANGLE- ' , 

F5.  1,  '  (M)  '  ) 

S  ( 1 , 1 ) =CG 
S(1,3)=-SG 
S (3,  1)=SG 
S  f  3 , 3) -CG 
ELSE 


WRITE  (LP.20)  COE .GAMMA 

20  FORMAT  (1H1 , /// , IX, ' C= ' , F4 . 2 , ' * (RP+RG) ;  INTERSECTING  ANGLE" ' , 
+  F5 . 1 ,  '  (M)  ' ) 

S (2 , 2) =CG 
S(2,3)=-SG 
S (3 , 2) =SC 
S (3 , 3) =CG 
END  IF 

C  SIMULATING  MESHING  OF  PINION  AND  GEAR (L=l  FOR  FINDING  INITIAL 
C  ROTATION  CORRESPONDING  TO  0  PINION  ROTATION)  AND  OUTPUT  RESULTS 
DO  25  L=1 , 2 
DO  35  1=1,4 
35  X(I)=0.D0 
DO  45  1=1, N 

X  (7) =FEEI*FLOAT (I- (N+l) /2) 

IF  (L.EQ.l)  X (7) =0. DO 
X  (5)  =  DARS  IN  (RP:tX  (7)  *SK/R) 

X (2) =X (7) /RMPG 
SF=DS IN (X (7) ) 

CF=DCOS (X (7) ) 

CALL  NONLIN 
X  (8)  =X  (2)  +X  (3) 

IF  (L.EQ.l)  THEN 
XIN=X (8) 

WRITE  (LP , 30) 

30  FORMAT  (////8X, 'FEE1 (D) ' ,8X, ' FEE2 (D) ' , 8X , ' K- ERROR (S) ' , 5X , 

'RP' ,13X, 'RG'/) 

GO  TO  25 
END  IF 

X (8) =X (8) -XIN 
X (7) =X (7) /DR 
X  (8)  =X  (8)  /DR 

X  (9)  =  (,X  (8)  -X  (7)  /RMPG)  *3600 .  DO 
Z  (I)  =X (8) 

ERR  (I)  =X  (9) 

WRITE  ( LP , 40)  (X(J) , J-7,11) 

40  FORMAT  (1X.5F15.7.F15.7) 

45  CONTINUE 

WRITE  (LP , 50) 

50  FORMAT  (//,’  FIND  THE  WORKING  RANGE  FOR  ONE  TOOTH:'/! 

DO  55  1=1, NT 

X  (7!  =  FEE  I ’-FLOAT  (I  -  (N*l)  /2)  /DR/2.D0 

X (8) =X  (7) *FII 

KK-I+NCOEF 


ANGLE=Z(KK)-Z(I) 

ERROR  (  I)  “(ANGLE-FII)  ’'3600. DO 
WRITE  ( LP , 60 )  X (7) ,X(8) .ANGLE 

60  FORMAT  (IX,  '  (’  ,F7.2,  * - '  ,F7.2,  ')  :  '  .F15.7.F15.7) 

55  CONTINUE  V-\ j 

do  65  1=1, nt 

ATEMP2=ERROR(I) 


IK? 


■  . 


K; 
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a  - 


;/ 


■* 
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GOTO  25 

75  IF  (DABS (ATEMP 1) .LT. DABS (ATEMP2))  1=1-1 
EMAX=0 . DO 
EMIN=0 . DO 
DO  85  J= 1 , NCOEF 
KS=I+J-1 
ET=ERR (KS) 

IF  (ET.LT.EKIN)  EMIN=ET 
IF  (ET.GT.EMAX)  EMAX=ET 
85  CONTINUE 

ET=EMAX~EMIN 

KK=I+NCOEF 

WRITE  (LP.80)  Z(I) ,Z(KK) ,ET 

80  FORMAT  (//IX, 'WORKING  RANGE  FOR  THE  GEAR  TOOTH:  ',F7.2,' - 

F7.2/1X, 'THE  MAXIMUM  KINEMATIC  ERROR:  ',F15.7,'  (S) 1 ) 

25  CONTINUE 
15  CONTINUE 

WRITE  (LP , 90) 

90  FORMAT ( 1 H 1 ) 

STOP 

END 

SUBROUTINE  FUNC 

THIS  SUBROUTINE  IS  USED  TO  GENERATE  FIVE  FUNCTIONS,  THAT  IS,  AT 
CONTACT  POINT,  POSITION  VECTOR  AND  NORMAL  OF  PINION  AND  GEAR 
MUST  COINCIDE 
IMPLICIT  REALi;8(A-H,0-Z) 

COMMON  /BLOCAl /  X ( 1 1 ) , Y(10) ,A(I0, 10) , Y1 (10) , IPVT(IO) .WORK (10) , 

EPS  I , DELTA, NC.NE.NDIM 

COMMON  / B LOCA2 /  ALPHA , RKS I , R P , RG , RMPC . R L . SA , CA , DR , CK , SK , HDC , C , 

S ( 3 , 3) , S  F , CF , SAK , CAK , Al , XP , YP , ZP , XG , YG , R , CONST (9) 
INPUT  PRE-DESIGN  ERROR  FUNCTION  EF  .AND  ITS  DIFFEVERTIVE  DEF, 
RECALLING  X(5)  IS  FEEP 
EF=-CONST ( 1) *X (5) *X (5) /2.D0 
DEF=-CONST(l)*X(5) 

CT=DCOS (X (A) ) 

ST=DSIN (X (A) ) 

CFEE=DC0S (X (3) ) 

SFEE=DS IN (X (3) ) 

CFEK=DCOS (X(3)*RKSI) 

SFEK=DS IN (X (3) *RKSI) 

GF=DARC0S  (  ( 1 .  DO+DEF*FMPG/  ( 1 .  DO+RMPG) )  *CK) 

CFGK=DCOS(X(5)*GF-RKSI) 

SFGK=DSIN(X(5)~GF-RKSI) 

CFP=DCOS (X (5) ) 

S  F  P  =  DS  I N  ( X  ( 5 ) ) 

CFPG=DCOS(X(5)-GF) 

SFPG=DS IN  ( X  (5)  *GF/' 

1 3  =  X i 5 ) /RMPG-EF^RKSI-GF 
XI  I  =  -  i  RP-RG)  '"SFP-RG':SFGK-RG-:T3:;CK"CFPG 
Yii-  !RP*RG;  !-CFP-RG,;CFGK»RG*T3,tCK*SFPG 
XP  =  X11*CT*R-'(CT-1.D0) 

Y  P  =  Y  1 1 

ZP=X  1 1  *ST*R*ST 
PNXP=CFPG’fCT 
?  NYP’SFPG 
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RNZP=CFDG*ST 

XG=RG*F  'EE+RG*X(2)*CK*CFEK 
YG=-RG*CFEE+RG*X  (2)  *CK*SFEK 
ZG=X  (1) 

RNXG-CFEK 

RNYG-SFEK 

RNZG-O.DO 

Y  (1)  =CF,VXP+SF*YP~S  (1,1)  *XG-S  (1 , 2)  *YG-S  (1 , 3)  *ZG 

Y (2) =CF*YP-SF*XP-S (2 , 1) *XG~S (2 , 2) *YG-S (2 , 3) *ZG~C 

Y  (3)  =ZP-S  (3,1)  *XG-S  (3 , 2)  *YG-S  (3 , 3)  *ZG 

Y  (4)  =CF,'CRNYP-SF*RNXP-S  (2 , 1)  *RNXG-S  (2 , 2)  *RNYG-S  (2 , 3)  *RNZG 
Y (5)  -RNZP-S (3 , 1) *RNXG-S (3 , 2) *RNYG-S (3 , 3) *RNZG 

X  ( 10)  =DSQRT  (XP*XP+ YP’^YP) 

X (1 1) =DSQRT (XG*XG+ YG*YG) 

C  WRITE  (6,20)  XP, YP,ZP,XG, YG,X(6) ,Al,A2,A3 
C  20  FORMAT  (IX, '$$$$’ .6F15.7) 

RETURN 

END 

C 

SUBROUTINE  NONLIN 
C 

C  THIS  SUBROUTINE  IS  USED  TO  SOLVE  NONLINEAR  EQUATIONS  BY  NEWTON- 
C  RAPHSON  METHOD 

IMPLICIT  REAL"8(A-H,0-Z) 

COMMON  /BLOCaI/  X (1 1) , Y (10) ,A(10,10) , Y1 CIO) , IPVT (10) ,WORK(10) , 
EPSI,DELTA,NC,NE,NDIM 

c 

DO  5  1=1, NC 
CALL  FUNC 

C  WRITE  (6,10)  I, (X(J) ,Y(J) , J-1.5) 

C  10  FORMAT  (IX,  '*>■*'  ,  15/5  (IX,  2D15 . 7/) ) 

DO  15  J= 1 ,NE 

IF  (DABS (Y(J)) .GT.EPSI)  GO  TO  25 
15  CONTINUE 
GO  TO  105 
25  DO  35  J-l.NE 
35  Y 1 ( J )  =  Y  ( J ) 

DO  45  J-l.NE 

DIFF=DABS  (X(J))  ’-DELTA 

IF  (X(J) .EQ.O.DO)  DIFF=DELTA 

XMAM=X ( J) 

x: J)=X(J)-DIFF 

Call  func 

X  J)=XMAM 
DO  55  K-l.NE 

55  A  >  K , J 1  - ( Y 1 (K)-Y(K)  )/DIFF 
45  CONTINUE 

DO  65  J-l.NE 
65  Y ' J ) = - Y 1  (J) 

CALL  DECOMP  (NDIM.NE.  A, CON’D,  IPVT,  WORK) 


SUBROUTINE  DECCHP  (NDIM , N , A, COND, IPVT , WORK) 

C 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  A  (NDIM, N) ,WORK(N) , IPVT(N) 

C 

C  DECOMPOSES  AREAL  MATRIX  BY  GAUSSIAN  ELIMINATION, 

C  AND  ESTIMATES  THE  CONDITION  OF  THE  MATRIX. 

C 

C  -COMPUTER  METHODS  FOR  MATHEMATICAL  COMPUTATIONS-,  BY  G.  E.  FORSYTHE, 
C  M.  A.  MALCOLM,  AND  C.  B.  MOLER  (PRENTICE-HALL,  1977) 

C 

C  USE  SUBROUTINE  SOLVE  TO  COMPUTE  SOLUTIONS  TO  LINEAR  SYSTEM. 

C 

C  IN  ll'T .  . 

C 

C  NDIM  =  DECLARED  ROW  DIMENSION  OF  THE  ARRAY  CONTAINING  A 

C  N  ORDER  OF  THE  MATRIX 

C  A  MATRIX  TO  BE  TRIANGULARIZED 

C 

C  OUTPUT.. 

C 

C  A  CONTAINS  AN  UPPER  TRIANGULAR  MATRIX  U  AND  A  PREMUTED 

C  VERSION  OF  A  LOWER  TRIANGULAR  MATRIX  I-L  SO  THAT 

C  (  PERML'l  ATION  MATRIX)  *A-L*U 

C 

C  COND  =  AN  ESTIMATE  OF  THE  CONDITION  OF  A. 

C  FOR  THE  LINEAR  SYSTEM  A*X  =  B  ,  CHANGES  IN  A  AND  B 

C  MAY  CAUSE  CHANGES  COND  TIMES  AS  LARGE  IN  X. 

C  IF  COND* 1.0  .EQ.  COND  ,  A  IS  SINGULAR  TO  WORKING 

C  PRECISION.  COND  IS  SET  TO  1.0D+32  IF  EXACT 

C  SINGULARITY  IS  DETECTED. 

C 

C  IPVT  *  THE  PIVOT  VECTOR 

C  IPVT(K)  =  THE  INDEX  OF  THE  K-TH  PIVOT  ROW 

C  IPVT(N)  •  (-1)** (NUMBER  OF  INTERCHANGES) 

C 

C  WORK  SPACE..  THU  VECTOR  WORK  MUST  BE  DECLARED  AND  INCLUDED 
C  IN  THE  CALL.  ITS  IN1UT  CONTENTS  ARE  IGNORED. 

C  ITS  OUTPUT  CONTENTS  ARE  USUALLY  UNIMPORTANT. 

C 

C  THE  DETERMINANT  OF  A  CAN  BE  OBTAINED  ON  OUTPUT  BY 
C  DET  f A)  =  IPVT (N  >  *  A(i.l)  *  A(2,2)  *  ...  *  A ( N , N )  . 

C 

IPVT rN )=1 

IF  ' N . EQ .  1 )  GO  TO  150 
NM1=N- 1 

C  COMPUTE  THE  1-NORM  OF  A  . 

ANORM=0.D0 
DO  20  J  =  1 , N 
T = 0 . DO 
DO  10  I  =  1 ,  N 
10  T=T-DABS (A f I , Jj  ) 

IF  fT.GT.ANOR.Ml  ANOR.M=T 
20  CONTINUE 

C  DO  GAUSSIAN  ELIMINATION  WITH  PARTIAL 

PIVOTING. 

DO  70  K= 1 , NM1 


?«»  f  «■ 


KP1-K+1 

C  FIND  THE  PIVOT. 

M=K 

DO  30  I-KPl.N 

IF  (DABS  (A  (I , K) )  .  GT .  DABS  (A (M ,  K)  )  )  M=I 
30  CONTINUE 
IPVT (K) =M 

IF  (M.NE.K)  IPVT (N) =-IPVT (N) 

t=a(m,k) 

a(m,k)=a(k,k) 

A  (K  ,  K)  =T 

C  SAIP  THE  ELIMINATION  STEP  IF  PIVOT  IS  ZERO 

IF  (T.EQ.O.DO)  GO  TO  70 
C 

C  COMPUTE  THE  MULTIPLIERS. 

DO  AO  I-KPl.N 
AO  A(I,K)=-A(I,K)/T 

C  INTERCHANGE  AND  ELIMINATE  BY  COLUMNS. 

DO  60  J=KP 1 , N 
T=A(M, J) 
a(m, j) =a(k, j) 
a(k, j)=t 

IF  (T.EQ.O.DO)  GO  TO  60 
DC  50  I-KPl.N 

50  A(I,J)-A(I,J)+A(I,K)*T 

60  CONTINUE 
70  CONTINUE 
C 

C  COND  »  (1-NORM  OF  A) * (AN  ESTIMATE  OF  THE  1-NORM  OF  A-INVERSE) 

C  THE  ESTIMATE  IS  OBTAINED  BY  ONE  STEP  OF  INVERSE  ITERATION  FOR  THE 
C  SMALL  SINGULAR  VECTOR.  THIS  INVOLVES  SOLVING  TWO  SYSTEMS 
C  OF  EQUATIONS,  (A-TRANSPOSE) *Y  =  E  AND  A"Z  =  Y  WHERE  E 
C  IS  A  VECTOR  OF  +1  OR  -1  COMPONENTS  CHOSEN  TO  CAUSS  GROWTH  IN  Y. 

C  ESTIMATE  =  (1-NORM  OF  Z) / (1-NORM  OF  Y) 

C 

C  SOLVE  (A-TRANSPOSE) *Y  =  E  . 

DO  100  K-l.N 
T=0. DO 

IF  (K.EQ.  1)  GO  TO  90 
KM1=K-1 
DO  80  1=1, KMl 
80  T  =  T  +  A(I,K)  ,vWORK  ( I ) 

90  EK= 1 . DO 

IF  (T.LT.O.DO)  EK=- 1 . DO 
IF  (A(K,K) .EQ.O.DO)  GO  TO  160 
100  WORK(K)  — (EK+T)/A(K,K) 

DO  120  KB= 1 , NM1 
K=N-KB 
T=0. DO 
KP 1 =K- 1 

DO  110  I-KPl.N 
110  T-T-A ( I , K  > -WORK (K) 

WORK ( K) =T 
M= I  PVT (K) 

IF  (M.EQ.K)  GO  TO  120 
T-WORK(M) 

WORK (M) -WORK (K) 

WORK(K)=T 


120  CONTINUE 
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C 

C 
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c 

c 
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c 
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c 
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c 

c 

c 
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YNORM=0 . DO 
DO  130  1=1, N 

130  YNORM=YNORM+DABS(WORK(I)) 

SOLVE  A*Z  =  Y 

CALL  SOLVE  (NDIM , N , A, WORK, IPVT) 

ZNORM=0 . DO 
DO  140  1=1, N 

140  ZNORM=ZNORM=DABS (WORK ( I) ) 


150 


ESTIMATE  THE  CONDITION. 

COND= ANORM*ZNORM/ YNORM 
IF  (CON'D.  LT.  1 .  DO)  C OND=1.DO 
RETURN 

1-BY-l  CASE. . 

COND= 1 . DO 

IF  (Ad,  Ij  .NE.O.DO)  RETURN 


160  COND=l .0032 
RETURN 


EXACT  SINGULARITY 


END 

SUBROUTINE  SOLVE  (NDIM. N . A, B , IPVT) 


IMPLICIT  REAL"8 ( A-H , 0~Z) 
DIMENSION  A (NDIM, N) , B (N) , IPVT(N) 


SOLVES  A  LINEAR  SYSTEM,  A*X  =  B 

DO  NOT  SOLVE  THE  SYSTEM  IF  DECOMP  HAS  DETECTED  SINGULARITY. 


-COMPUTER  METHODS  FOR  MATHEMATICAL  COMPUTATIONS-,  BY  G.  E.  FORSYTHE, 
M.  A.  MALCOLM,  AND  C.  B.  MOLER  (PRENTICE-HALL,  1977) 

IN1UT.  . 

NDIM  =  DECLARED  ROW  DIMENSION  OF  ARRAY  CONTAINING  A 
N  =  ORDER  OF  MATRIX 

A  =  TR I ANGULAR IZED  MATRIX  OBTAINED  FROM  SUBROUTINE  DECOMP 

B  =  RIGHT  HAND  SIDE  VECTOR 

IPVT  =  PIVOT  VECTOR  OBTAINED  FROM  DECOMP 

OUTPUT. . 


B  =  SOLUTION  VECTOR,  X 


DO  THE  FORWARD  ELIMINATION. 

IF  ON. EQ. 1 )  GO  TO  50 
NM1=N-  1 
DC  20  K= 1 , NM 1 
KPi=K=l 
M= I  PVT (K) 

T=B (M) 

B(M)=B(K) 

B (Kj  =T 

DO  10  I =KP 1 , N 


SB 

ij 

I 


a 


i 

i 


72 


£ 


10  B(I)=B(I)+A(I,K)*T 

20  CONTINUE 

NOW  DO  THE  BACA  SUBSTITUTION 

DO  40  KB=1,NM1 
KM1=N-KB 
K=KM1+1 

b(k)=b(k)/a(k,k) 

T— B(K) 

DO  30  1=1, KM1 
30  B(I)=B(I)+A(I,K)*T 

40  CONTINUE 
50  B(1)=B(1)/A(1,1) 

RETURN 

END 


FLOWCHART  FOR  PROGRAM  II 


c. . 
c. . 
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PROGRAM  I I  * 

SURFACE  OF  PINION  GENERRATED  BY  CONE  CUTTER  * 

* 

AUTHORS:  FAYDOR  LITVIN  * 

JIAO  ZHANG  * 


*  *  *  *  *  *  it  it  it  it  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  i(  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  *  it  it  it  it  it  it  it  it  it  it  it  it 


PURPOSE 

THIS  PROGRAM  IS  USED  TO  CALCULATE  THE  SURFACE  OF  A  PINION  WHICH 
GENERATED  BY  CONE  CUTTER 


HIS  PROGRAM  IS  WRITTEN  IN  FORTRAN  77.  IT  CAN  BE  COMPILED  BY  V 
COMPILER  IN'  IBM  MAINFRAME  OR  FORTRAN  COMPILER  IN  VAX  SYSTEM. 


c , 

C. 

C. 

c. 
c. 
c. 
c. 
c. 
c. 
c 

r 

c 
c 
c 
c 

c  NOTE 
C 
C 
C 
C 


IMPLICIT  REAL”8 (A-H.O-Z) 

DIMENSION  NS (15) ,XP(15,999) ,YP(15,999) ,XERR0R(999) 

C 

C  DEFINE  PARAMETERS  USED  BY  PROGRAMS 

C 

C  (l:  IN  ANG  LP  .ARE  UNIT  NUMBERS  ASSIGNED  TO  THE  INPUT  AND  OUTPUT 
C  DEVICES 


IN-5 


IS 


LP=6 

C  (2)  NDBUG  IS  USE  TO  CONTROLS  THE  AUXILIARY  OUTPUT  FOR  DEBUGGING 
NDBUG-2 

C  <'3)  NC  IS  THE  UPPER  LIMITATION  OF  REPEATATION  FOR  SOLVING  NONLINEAR 
C  EQL'TIONS ; 

C  EPSI  IS  THE  CLEARANCE  OF  FUNCTION  VALUES  WHEN  THE  FUNCTIONS 

C  IS  CONSIDERED  AS  SOLVED  (ALL  FUNTIONS  HAVE  FORMS  OF  F(X)-O); 

C  RKSTA  IS  THE  RELATIVE  DIFFERENCE  FOR  TAKING  DERIVATIVES 

C  NC,  EPSI  AND  RKSTA  MAY  BE  CHANGED  WHEN  SOLUTIONS  ARE  DIVERGENT 

C  OR  LESS  ACCURATE 

NC= 1 00 
RKSTA-l.D-3 
EPS  I -  1 . D- 1 2 

C  Li  OTHER  PARAMETERS  (DON"  T  CHANGE; 

DP  =  DAT  AN  ( I . DO ) / L  5 . DO 


r 

C  DEFINE  INPUT  PAR AMTERS  OF  PROBLEM (USE  INCH  AS  UNIT  OF  LENGTH) 


C  I  PINION  AND  GEAR:  PN  =  DI AMETRAL  PITCH:  Nl-PINiON  TOOTH  NUMBER; 

C  RilS-PRESSURE  ANGLE  '.DEGREE;  ;  HDC-HEIGHT  OF  DEDENDUM  OF  PINION; 

G  HAC-HEIRHT  OF  ADDENDUM  OF  PINION 

PN-10.D0 
N  1  -  2  0 

R KS  =  20 .  D0';'DR 
HDC= I . DO/PN 
HAC  I .DO/PN 

2  TOOL:  ALPHA-HALF  OF  CONE  VERTEX  ANGLE ( DEGREE) ; 

PC-RADIL'S  OF  BOTTOM  CIRCLE  OF  CONE 


V.  -vv.-,  V  .■  'j-'/jV,  no. 

>  J ►  A  .  •  /•  -V  .  *  JV  h  ■  *■  _N-  A 


ALP=80.D0*DR 
RC-1 .DO 

(3)  MISALIGNMENT:  NMIS=ID  NO. (l=CROSSING  AXES,  2= INTERSECTING  AXES); 
NG'NO.  OF  MISALIGNED  ANGLES  TO  BE  SIMULATED  (FROM  0.  TO 
(NG-1)*GAMMAI) ;  GAMMAI= INCREMENT  OF  MISALIGNED  ANGLE (MINUTE) ; 
N=500 
NG=2 

GAMMAI  =  5  .DO 

(A)  OUTPUT:  ZI“INCREMENT  OF  TOOTH  LENTGH  OF  CROSS  SECTION; 

NL=NO.  OF  CROSS  SECTIONS  WHERE  PINION  PROFILE  IS  .sinuLATED; 

N=NO .  OF  POINTS  USED  TO  DEVIDE  TOOL  PROFILE  IN  ANY  CROSS 
SECTION (LARGE  N,  MORE  POINTS  IS  GOT  FOR  PINION  TOOTH  PROFILE) 
ZI=0. 10D0 
NL=3 
N=500 


C  DESCRIPTION  OF  OUTPUT  PARAMETERS 
C 

C  Z1=DISTANCE  BETWEEN  CROSS  SECTION  CONSIDERED  AND  MIDDLE  CROSS 

C  SECTION 

C  NO=OUTPUT  NO. 

C  Y i=TOOL  SURFACE  AUXILIARY  VARIABLE 

C  XP=X  COORDINATE  OF  PINION  PROFILE 

C  YP= Y  COORDINATE  OF  PINION  PROFILE 

C  R1=RADIUS  OF  PINION  PROFILE 

C  FEE=CORESPONDING  PINION  SURFACE  PARAMETER 

C  XSH=AVERAGE  DEVIATION  SHIFT  OF  CROSS  SECTION  PROFILE  FROM  PROFILE 

C  OF  MIDDLE  SECTION (INVOLUTE  CURVE) 

C  XPE-MAXIMUM  DEVIATION  OF  CROSS  SECTION  PROFILE  FROM  INVOLUTE  CURVE 

C  SDX= STANDARD  DEVIATION  OF  CROSS  SECTION  PROFILE  FROM  INVOLUTE 

C  CURVE 


FIND  AUXILIARY  VALUES  FOR  CALCULATION 
RP=FL0AT(N1) /2.D0/PN 
CL=RC/DSIN (ALP) 

D=CL’;DC0S  (ALP) 

Al=CL-HDC/DCOS(RKS) 

RPU=RP*HAC 
DO  5  1=1 , N  L 
Z  =  FLOAT  ( I - 1 )  *Z  I 
YY=D-Z/DTAN (ALP) 

WRITE  (LP.10)  Z 

10  FORMAT  (1H1/1X, 'Zl='  .F15.7/1X,  'NO' ,10X. ’ Y 1 '  .  13X,  ' XP ' , 13X,  ' YP '  , 
13X, 'R1 ' , 13X, ' FEE' ) 

CALCULATE  THE  PROPILE  OF  THE  SURFACE  CUT  BY  PLANE  Z=CONST 
KKK=0 

DO  15  J  =  1 .  N 

Y I  =  Y  Y ’T LOAT <  J -  1 : /FLOAT  ! N- 1 ) 

F-DSOFT • l.-'Z ' (D-Yl ) /DTAN (ALP) ) **2) 

"EE=  'D-Yl  ■  ”F/DC0S  '  ALP  )  -  A  i F“DCOS  (ALP )  "DCCS  t  ALP  )  *DS  IN  ( ALP  1  *DSIN  ( 
ALP  ■ -RP/ (DSIN1ALP • *DCOS (ALP-RKS) -F*DCOS (ALP) ”DSIN (ALP-RKS) ■ 
XP  1 , :  = : D-Yl ) " (DTAN (ALP)  *‘F*DCOS (ALP-RKS- FEE J ~DS IN (ALP-RKS ‘FEE )  > 

- Al "DS IN ( FEE-RKS )  -RP’!FEE"DCOS  i  FEE)  ‘RP*DS IN  (  FEE ) 
vp  i  ,  ■  ...  i D-Yl )  ’•  (DTAN  (ALP)  *F*DSI.N  (ALP-RKS* FEE)  *DCOS( ALP-RK.S* FEE) ) 
-A  I  DC  OS  (  FEE-RKS)  -RP:rFEE:'DS  IN' '  FEE  !  *RP’-DC0S  <  FEE) 

?  1’DSQR T  ■  XP  (  I  .  J  )  *  +  2~  YP  ( I .  J )  :r ’■ 2  ) 

IF  ' F  1 . GT . RPU )  THEN 


N.V.V..  ■ 

„y.  o-  >.v; 


v.v, 


XP(I,  J)=XP(I,  JJ)  +  (XP(I,  J)-XP(I,  JJ))/(R1-R1TEMP)’,'(RPU-R1TEMP) 

YP (I , J)=YP (I , JJ) + (YP (I , J)-YP (I , JJ)) / (R1-R1TEMP) * (RPU-R1TEMP) 
FEE=FEETEM+ (FEE-FEETEM) / (R1-R1TEMP) * (RPU-R1TEMP) 

Y 1= Y1TEM* ( Yl- Y ITEM) / (R 1-R 1TEMP) * (RPU-R 1TEMP) 

R1=DSQRT (XP (I , J) **2+YP (I , J) **2) 

KKK=1 
END  IF 

WRITE  (LP.20)  J,Y1,XP(I, J) ,YP(I, J) ,R1,FEE 
20  FORMAT  (IX, I4.4F15.7.F15.7) 

IF  (KKK.CT.O)  GO  TO  30 
FEETEM=FE£ 

Y1TEM=Y1 
R1TEMP=R1 
15  CONTINUE 
30  NS (I) = J- 1 

IF  (I.NE.  1)  GO  TO  55 
NSl-NS(l) 

GO  TO  5 

C  PREPARATION  OF  INTERPLORATION 
55  NS2=NS ( I ) 

DO  105  L-1.NS2 
J=NS2* i-L 

IF  (YPd.NSl)  .GT.YP(I,  J)1  GO  TO  110 
105  CONTINUE 
110  NS2  =  J 
NREC=2 
XERS=0 . DO 
DO  115  L-1.NS2 
DO  125  J=NREC , NS  1 
IF  (YP(1, J) .GT.YP(I.L) )  GO  TO  120 
125  CONTINUE 
120  NREC=  J 
Jl  =  J- 1 

XERROR  (LI  =  XP  (I  ,  L)  ~XP  (1 ,  J 1 )  -  (YP  (I ,  L)  -  YP  (1 ,  Jl) )  *  (XP  (1 ,  Jl  -XP  (1 ,  J1 )  1 
/(YP(1, J)-YP(1,  Jl)) 

XERS=XERS*XERROR (L) 

115  CONTINUE 

XERS=-XERS /FLOAT (NS2) 

XPE=0 . DO 
SDX=0 . DO 

IF  (NDBUG.GT.2)  WRITE  (LP,80) 

80  FORMAT  (1H1.13X. 'NO. ' ,4X. 'DIVIATION  VALUE') 

DO  45  L=1  ,NS2 

XERROR (L) =XERROR (L)-XERS 

IF  (NDBUG.GT.2)  WRITE  (LP.50)  L , XERROR ( L ) 

50  FORMAT  (13X. I3.2F15 .7) 

IF  (DABS (XERROR (L) ) .GT.XPE;  XPE=DABS (XERROR ( L ) ) 

SDX=SDX-XERROR (L) **2 
45  CONT’NL'E 

SDX=DS0RT(SDX/FL0AT(NS2) ) 

WRITE  (LP.40)  XERS.XPE.SDX 

40  FORMA.T  '///IX.  ' XSH= ' .E15.7.5X, 'XPE=' .E15.7.5X,  ' SDX= ' . El 5. 7) 

IF  'NDBUG.GT.2)  WRITE  (LP.60)  NS1.NS2 


60  FORMAT  f //5X. ’NS1= ' , 16, 6X, 'NS2= ’ , 16) 
5  CONTINUE 
STOP 
END 
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PROGRAM  III 

UDDERCUTTING  CONDITION  FOR  PINION  GENERRATED 
BY  CONE  CUTTER 


AUTHOR:  FAYDOR  LITVIN 
JIAO  ZHANG 


it  it  it  it  it  it  it  it  it  it  A  it  it  it  it  it  it  it  A  A  A  A  A  A  it  it  A  it  it  it  it  it  it  it  it  it  A  A  A  A  A  A  A  A  it  A  A  A  A  A  A  A  A  A  A 


PURPOSE 

THIS  PROGRAM  IS  USED  TO  FIND  THE  UDDERCUTTING  CONDITIONS  FOR  A 
PINION  GENERATED  BY  CONE  CUTTER 


THIS  PROGRAM  IS  WRITTEN  IN  FORTRAN  77.  IT  CAN  BE  COMPILED  BY  V 
COMPILER  IN  IBM  MAINFRAME  OR  FORTRAN  COMPILER  IN  VAX  SYSTEM. 

IMPLICIT  REAL’-8(A-H,0-Z) 

DEFINE  PARAMETERS  USED  BY  PROGRAMS 

(1)  IN  ANG  LP  ARE  UNIT  NUMBERS  ASSIGNED  TO  THE  INPUT  AND  OUTPUT 

DEVICES 

IN=5 

LP=6 

(2)  NDBUG  IS  USE  TO  CONTROLL  THE  AUXILIARY  OUTPUT  FOR  DEBUGGING 
NDBUG=2 

(3)  OTHER  PARAMETERS (DON ' T  CHANGE) 

DR=DATAN ( 1 . DO) / 45 . DO 

DEFINE  INPUT  PARAMTERS  OF  PROBLEM (USE  INCH  AS  UNIT  OF  LENGTH) 

(1)  PINION  AND  GEAR:  PN=DI AMETRAL  PITCH;  Nl=PINION  TOOTH  NUMBER; 

RKS I =PRESSURE  ANGLE;  HDC=HEIGHT  OF  DEDENDUM  OF  PINION; 

TL=TOOTH  LENGTH 
PN= 10 . DO 
Nl  =  20 

RKSI=20.D0*DR 
HDC  = 1 .DO/PN 
TL=6 . DO/PN 

(2)  TOOL:  ALPHA=HALF  OF  CONE  VERTEX  ANCLE (DEGREE) ; 

RC=RADIUS  OF  BOTTOM  CIRCLE  OF  CONE 
ALPHA=80.  DO’-DR 
RC=1 .DO 

(3)  PROBLEM:  NPROB=ID  NO.  OF  PROBLEM  (-1=GIVEN  N1  AND  HDC ,  FIND  IF 

UDDERCUTTING  OCCUR;  0=GIVEN  Nl,  FIND  MAXIMUM  HDC  WITHOUT  UDDER¬ 
CUTTING;  1=GI VEN  HDC,  FIND  MINIMUM  Nl  WITHOUT  UDDERCUTTING); 
N=NO .  OF  THETA  VALUES  USED  CALCULATION  (BETWEEN  THETAS  CORRES¬ 
PONDING  TO  MIDDLE  SECTION  AND  EDGE  SECIION  OF  PINION) 

NPROB=  1 
N=  1 1 
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C  DESCRIPTION  OF  OUTPUT 
C 

C  OUTPUT  IS  A  STATEMENT  BASED  ON  THE  PROBLEM  WITHOUT  ANY  LITERAL 
C  PARAMETER 
C 

C  FIND  AUXILIARY  VALUES  FOR  CALCULATION 
CL=RC/DSIN (ALPHA) 

SA=DS IN (ALPHA) 

CA=DCOS (ALPHA) 

SK=DS IN (RKS I) 

CK=DC0S (RKSI) 

SAK=DS IN ( ALPHA-RKS I ) 
cak=dcos (alpha-rksi) 

THEMAX=DARSIN(TL/2.0/SA) 

IF  (NPROB.NE.O)  THEN 
Al= (CL-HDC/CK) *CK 
ACL=Al/CL 
END  IF 

IF  (NPROB.LT. 1)  THEN 
RP= FLOAT (Nl)/2. /PN 
RCL=RP/CL 
END  IF 

WRITE  (LP.90) 

90  FORMAT  ( 1H1 ) 

IF  (NPROB)  5,15,25 
C  CHECK  IF  UNDERCUTTING  OCCURS 

5  IF  (NDBUG.LT. 1)  WRITE  (LP,10) 

10  FORMAT  (1H1/3X, 'NG' ,7X, ' THETA' , 12X , ' A ' , 13X, ' B ' , 1  AX, ' C ' , 10X , 

'  B*>:2-4.  *A*C  '  ,3X,  'U1/(RC/S  IN  (ALPHA!)  ') 

UMIN=5 . DO 
DO  45  1=1 , N 

THE=THEMAX*DBLE(FLOAT(I-l) ) /DBLE (FLOAT (N- 1) ) 

ST=DSIN (THE) 

CT=DCOS (THE) 

C0E1=CA*CAK+CT*SA*SAK 

COE2=SA*CAK-CT*CT*CT*CA*SAK 

COE3=SA*CAK-CT*CA"SAK 

COE4=ST,!ST*SA*CA 

A=COEl*COE2 

B=-aCL,v(COE2+COE1*COE4)-RCL::COE3-:*3 

C=ACL"ACL*C0E4 

DD=B*B-4.*A*C 

IF  (DD.GE.O.)  THEN 

UlCL=(-B-DSQRT(DD))/2./A 

IF  (UMIN.GT.U1CL)  UMIN=U1CL 

t’2CL=  (-B-DSQRT  (DD) )  12.1  h 

IF  (NDBUG.LT. 1)  WRITE  (LP.100)  I , THE. A, B , C . DD. U1CL 
100  FORMAT  (IX, I4.8F15.7) 

ELSE 

IF  (NDBUG.  L'l  .  1)  WRITE  (LP.100)  I .  THE .  A .  B  ,  C  ,  DD 
END  IF 
45  COMINUE 

IF  (UMIN.LT. 1.D0)  WRITE  (LP , 110) 

110  FORMAT  (///IX, 'UDDERCUTTING  WILL  OCCUR  FOR  YOUR  DESIGN’) 

IF  (UMIN.GE. l.DO)  WRITE  (LP.120) 

120  FORMAT  (///IX, 'UDDERCUTTING  WILL  NOT  OCCUR  FOR  YOUR  DESIGN') 
GO  TO  35 

C  DETERMINE  THE  MAXIMUM  ADDENDUM  HEIGHT  OF  RACA  CUTTER 
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15  IF  (NDBUG.LT. 1)  WRITE  (LP.20) 

20  FORMAT  (1H1/3X, 'NO' ,7X, 'THETA' , 12X, 'A' , 13X, 'B' , 14X, 'C' , I OX, 

+  ’B**2-4.*A*C' ,3X, 'ALLOWED  RATIO  OF  HDC/(1/PN)') 

UMIN=5 . DO 
DO  55  1=1, N 

THE=THEMAX*DBLE (FLOAT (1-1) ) /DBLE (FLOAT (N-l) ) 

ST=DSIN (THE) 

CT=DCOS (THE) 

C0E1=CA*CAK+CT*SA*SAK 

COE2-SA*CAK-CT*CT*CT*CA*SAK 

COE3=SA*CAK-CT*CA*SAK 

C0E4=ST*ST*SA*CA 

C=COE 1 *COE2-RCL*COE3**3 

B=-(C0E2+C0E1*C0E4) 

A=C0E4 

DD=B*B-4.  *A*C 
IF  (DD.GE.O.)  THEN 
IF  (A.NE.O.)  THEN 
U1CL- (-B-DSQRT (DD) ) /2 . /A 
U2CL- (-B-DSQRT (DD) ) /2 . /A 
ELSE 

U1CL=-C/B 

U2CL-0. 

END  IF 

U1CL-(CK-U1CL)*CL*PN 
L'2CL=  (CK-U2CL)  *CL*PN 

IF  (NDBUG.LT. 1)  WRITE  (LP.100)  I , THE , A, B , C , DD , U 1CL 
ELSE 

IF  (NDBUG.LT. 1)  WRITE  (LP.100)  I , THE , A , B , C , DD 
END  IF 

IF  (UMIN.GT.U1CL)  UMIN-U1CL 
55  CONTINUE 

WRITE  (LP.200)  UMIN 

200  FORMAT  (///IX, 'TO  AVOID  UDDERCUTTING,  IT  IS  NECESSARY  TO  KEEP  DEDE 
-NDUM  OF  PINION  <= ' , F10 . 7 , ' /PN ' ) 

GO  TO  35 

C  DETERMINE  THE  MINIMUM  NO.  OF  TEETH  FOR  UNUNDERCUTTING 
25  IF  (NDBUG.LT. 1)  WRITE  (LP,30) 

30  FORMAT  (1H1./3X,  'NO'  ,7X.  'THETA'  ,7X,  'NO.  OF  TEETH') 

RNMAX-0. 

DO  65  1=1 , N 

THE  =  THEMAX,rDBLE (FLOAT ( I —  1 ) ) /DBLE (FLuAT (N-l) ) 

ST=DS IN (THE) 

CT=DCOS (THE) 

C0E1=CA*CAK+CT*SA*SAK 
COE2-SA*CAK-CT*CT*CT*CA*SAK 
CCE3=SA*CAK-CT*CA*SAK 
C  0  E  4  =  S  T ::  S  T  *  S  A C  A 

RR=(C0E1*C0E2-ACL"(C0E2*C0E1,';C0F4)-ACL"ACL"C0E4)  /COE3**3 
R N  =  2  .  ■•RR*PN*CL 
IF  (RNMAX.LT. RN)  rnmax=rn 
IF  (NDBUG.LT. 1)  WRITE  (LP.100)  I.THE.RN 
65  CONTINUE 

WRITE  (LP , 300)  RNMAX 

300  FORMAT  (///  IX, ’WITHOUT  UDDERCUTTING,  MINIMUM  TOOTH  NO.  OF  PINION 
-IS: ', FI  1.7) 

35  WRITE  ( LP , 400) 

400  FORMAT ( 1 H 1 ) 
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C  PURPOSE 
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PROGRAM  IV 

CONTACT  ELLIPSIS  FOR  PINION  GENERATED  BY  CONE 
CUTTER  IN  MESHING  WITH  REGULAR  GEAR 

AUTHORS:  FAYDOR  LITVIN 
JIAO  ZHANG 


Ve  it  *  :Y  *  *  Vr  *  *  *  *  *  *  *  it  it  *  *  *  *  it  *  Vr  *  *  *  *  Vt  *  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it 


THIS  PROGRAM  IS  USED  TO  CALCULATE  THE  CONTACT  ELLIPSIS  WHEN  PINION 
GENERATED  BY  CONE  CUTTER  IN  MESHING  WITH  REGULAR  INVOLUTE  GEAR 


P.-.' 


m 


c 

C  THIS  PROGRAM  IS  WRITTEN  IN  FORTRAN  77.  IT  CAN  BE  COMPILED  BY  V 
C  COMPILER  IN  IBM  MAINFRAME  OR  FORTRAN'  COMPILER  IN  VAX  SYSTEM. 

C 

IMPLICIT  REAL*8 (A-H,0-Z) 

C 

C  DEFINE  PARAMETERS  USED  BY  PROGRAMS 
C 

C  (1)  IN  ANG  LP  ARE  UNIT  NUMBERS  ASSIGNED  TO  THE  IN1UT  AND  OUTPUT 
C  DEVICES 

I  N= 5 
LP=6 

C  (2)  NDBUG  IS  USE  TO  CONTROLL  THE  AUXILIARY  OUTPUT  FOR  DEBUGGING 
NDBUG=2 

C  (3)  OTHER  PARAMETERS (DON'T  CHANGE) 

DR  =  DATA\'(1.D0)/45.D0 
C 

C  DEFINE  INPUT  PARAMTERS  OF  PROBLEM (USE  INCH  AS  UNIT  OF  LENGTH) 

C 

C  (1)  PINION  AND  GEAR:  PN=DIaMETRAL  PITCH;  Nl=PINION  TOOTH  NUMBER; 

C  RMPG=TOOTH  NUMBER  RAT  10 (GEAR  TOOTH  N0./N1); 

C  RKS I=PRESSURE  ANGLE 

PN= 10. DO 
N I  =  2  0 
RMPG=2 . DO 
RKS I  =  20 . DO"DR 

C  (2)  TOOL:  ALPHA=HALF  OF  CONE  VERTEX  ANGLE (DEGREE) ; 

C  RC-RADIUS  OF  BOTTOM  CIRCLE  OF  CONE; 

ALPHA"89.5D0'-'DR 
RC= 1 . DO 

C  (3)  DEFORMAT  I  ON :  DEL=CONTACT  DEFORMATION  AT  CONTACT  POINT 
DEL=A.D-4 

C  (4)  OUTPUT:  PD=REC I PROCAL  OF  INCREMENT  OF  S  VALl’EiS  IS  THE  DISTANCE 
C  CONTACT  POINT  AND  INSTANTANEOUS  CENTER) 

PD= 1 00. DO 
C 

C  DESCRIPTION  OF  OUTPUT  PARAMETER 

C 

C  S=DI STANCE  BETWEEN  CONTACT  POIN"  -.ND  INSTANTANEOUS  CENTER 
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C  Rl-PINION  RADIUS  OF  CONTACT  POINT 

C  FEE=PINION  ROTATION  ANGLE  CORRESPONDING  TO  CONTACT  POINT 
C  A= LENGTH  OF  HALF  SHORT  AXIS  OF  CONTACT  ELLIPSE 

C  B=LENTH  OF  HALF  LONG  AXIS  OF  CONTACT  ELLIPSE (ALONG  DIRECTION  OF 
C  GEAR  TOOTH  LEHGTH) 

C 

C  FIND  AUXILIARY  VALUES  FOR  CALCULATION 
RP= FLOAT (Nl)/2./PN 
SA=DS IN (ALPHA) 

CA^DCOS (ALPHA) 

SK=DSIN (RKSI) 

CK=DCOS (RKSI) 

RL=RC/SA 

SMAX= (DSQRT ( ( ( 1 . DO+2 . DO/ FLOAT (Nl) ) /CK) **2- 1 . DO) -SK/CK) *CK*RP 
NA=IDINT(SMAX*PD+0.5D0) 

SMAX=FLOAT(NA)/PD 

SMIN= (DSQRT (.(1.-2 . /FLOAT (Nl) /RMPGl /CK) **2~ 1 . ) -SK/CK) *CK*RP*RNPG 
NB=- IDINT (SMIN*PD+0 . 5D0) 

SMIN=FLOAT(NB) /PD 
N-l+NA-NB 

WRITE  (LP.IO)  RP , N 1  , RMPG 

10  FORMAT  (1H1 ,///5X, 'RADIUS  OF  THE  PITCH  CIRCLE  OF  PINION  (RP) ; ' , 
-F15. 7//5X, ' TOOTH  NO.  OF  PINION (Nl) I8//5X, ' RATIO  OF  ANGULAR  VELO 
■►CITY  (OMEGA  OF  PINION/OMEGA  OF  GEAR) F15 . 7///1X, ' OUTPUT  DESIGNAT 
-ION:  S=DI STANCE  BETWEEN  CONTACT  POINT  AND  INSTANTANEOUS  CENTER'/ 
-21X, 'Rl-PINION  RADIUS  OF  CONTACT  POINT;  FEE^PINION  ROTATION  CORRES 
-PONDING  TO  CONTACT  POINT ’ /21X, ' A=LENGTH  OF  HALF  SHORT  AXIS:  B=LENG 
-TH  OF  HALF  LONG  AXIS (ALONG  DIRECTION  OF  GEAR  TOOTH  LEHGTH)'///// 
-10X, 'S ' , 14X, 'R1 1 , 10X, ' FEE1 (D) 1 , 1 IX , ' A ' , 14X , ' B ' , 1 3X , ' B/A ' //) 

C  FIND  POSITION  AND  DIMENSION  OF  CONTACT  ELLIPSE 
DO  5  L=1,N 

S=SMIN- (SMAX-SMIN) /FLOAT (N-l)-FLOAT(L-l) 

A=DSQRT  (DEL*RP*  (SK-S/RP)  *  (RMPG*SK-S/RP)  *2  .  DO/'  ( 1 .  DO-RMPG)  /SK) 
B-DSQRT(DEL*2.D0*Sa/CA*  (RL-2.D0*RP/CK/FL0AT(N1)-S*SK/CK) ) 

C  =  B/A 

FEE-S/RP/CK 

R  1=RP*DSQRT  ( 1 .  DO-  (FEE':CK)  **2+2 .  DO*FEE*CK*SK) 

FEE=FEE/DR 

WRITE  (LP , 20)  S , R 1 , FEE . A, B , C 
20  FORMAT ( IX , 6F 1 5 . 7) 

5  CONTINUE 
WRITE  (LP , 30) 

30  FORMAT  ( 1 H 1 ) 

STOP 

END 
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*  PROGRAM  V 

*  KINEMATIC  ERROR  OF  THE  PINION  GENERATED  BY  CONE 

*  CUTTER  MESHING  WITH  A  MISALIGNED  REGULUAR  GEAR 


AUTHORS:  FAYDOR  LITVIN 
JIAO  ZHANG 


*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  V;  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 


PURPOSE 


THIS  PROGRAM  IS  USED  TO  CALCULATE  THE  KINEMATIC  ERROR  OF  A  PINION 
GENERATEED  BY  CONE  CUTTER  MESHING  WITH  A  MISALIGNED  REGULAR 
INVOLUTE  GEAR 


THIS  PROGRAM  IS  WRITTEN  IN  FORTRAN  77.  IT  CAN  BE  COMPILED  BY  V 
COMPILER  IN  IBM  MAINFRAME  OR  FORTRAN  COMPILER  IN  VAX  SYSTEM. 


IMPLICIT  REAL "8 (A_H ,  0~Z) 

DIMENSION  Z (99 ) , ERR (99) .ERROR (99) 

COMMON  /BLOCAl/  X ( 1 1) , Y ( 10) , A ( 10 , 10) , Y 1 (10) . IPVT ( 10) , WORK ( 10) 
EPSI .DELTA, NC.NE.NDIM 

CC  MON  /BLOCA2/  ALPHA, RKSI.RP.RG.RMPG.RL.SA.CA. DR. CK.SK.HDC, 

C , S (3 , 3) .SF.CF.SAK.CAK.Al.XP, YP.ZP.XG, YG.R 


DEFINE  PARAMETERS  USED  BY  PROGRAMS 


il)  IN  ANG  LP  ARE  UNIT  NUMBERS  ASSIGNED  TO  THE  INPUT  AND  OUTPUT 
DEVICES 

IN-5 
LP  =  6 

2)  NDBUG  IS  USE  TO  CONTROLL  THE  AUXILIARY  OUTPUT  FOR  DEBUGGING 
NDBUG-2 

(3)  NC  IS  THE  UPPER  LIMITATION  OF  REPEATATION  FOR  SOLVING  NONLINEAR 
EQUTIONS; 

EPSI  IS  THE  CLEARANCE  OF  FUNCTION  VALUES  WHEN  THE  FUNCTIONS 
IS  CONSIDERED  AS  SOLVED  (ALL  FUNTIONS  HAVE  FORMS  OF  F(X)-O); 
DELTA  IS- THE  RELATIVE  DIFFERENCE  FOR  TAKING  DERIVATIVES 
NC.  EPSI  AND  DELTA  MAY  BE  CHANGED  WHEN  SOLUTIONS  ARE  DIVERGENT 
OR  LESS  ACCURATE 
NC= 100 
DELTA- 1 . D-3 
EPS  1= I .D- 12 

l  other  parameters (don't  change) 

NDI.M-10 

NF.-5 

DR-DATANM  .DO'  /45.D0 

DEFINE  INPUT  PARAMTERS  OF  PROBLEM ( USE  INCH  AS  UNIT  OF  LENGTH) 

I!  PINION  AND  GEAR:  PN=DlAMETRAL  PITCH;  Nl-PINION  TOOTH  NUMBER; 
RMPG-TOOTH  NUMBER  RATIO (GEAR  TOOTH  NO./N1); 

RKS I -PRESSURE  ANGLE;  HDC=HEIGHT  OF  DEDENDUM  OF  PINION; 
COE-COEFF.  OF  CENTRAL  DISTANCE (USUALLY  COE= 1 . ) 


I*.' 


L* 


UNCLASSIFIED 


SPUR  SEARS'  OPTIHAL  GEONETRV  METHODS  ?0R TiENERATION 
TOOTH  CONTACT  ANA  (U)  ILLINOIS  UNIV  AT  CHICAGO  CIRCLE 
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Nl-20 
RMPG-2.D0 
RKS I  =20 .  D0’cDR 
HDC-1 .DO/PN 
C0E=1 . 005D0 

(2)  TOOL:  ALPHA-HALF  OF  CONE  VERTEX  ANGLE (DEGREE) ; 

RC-RADIUS  OF  BOTTOM  CIRCLE  OF  CONE 

R-RADIUS  OF  ARC (FOR  EXACT  CONE  CHOOSE  R  >1 .D6*Nl/2/PN) 

ALPHA-80. do*dr 
RC-1 .DO 
R-1.0D3 

(3)  MISALIGNMENT:  NMIS-ID  NO. (1-CROSSING  AXES,  2- INTERSECTING  AXES); 

NG-NO.  OF  MISALIGNED  ANGLES  TO  BE  SIMULATED  (FROM  0.  TO 
(NG-1) *GAMMAI) ;  GAMMA I- INCREMENT  OF  MISALIGNED  ANGLE (MINUTE) ; 
NMIS-1 
N'G=2 

GAMMAI-5  .  DO 

v4)  OUTPUT:  FEE  I  =  INCREMENT  OF  ROTATION  ANGLE  OF  PINION (DEGREEE) 

FEE  I- 1 .0D0:;DR 

DESCRIPTION  OF  OUTPUT  PARAMERTERS 

FEE1 -ROTaT I  ON  ANGLE  OF  PINION 
FEE2— ROTAT ION  ANGLE  OF  GEAR 
RP-RADIUS  OF  PINION  CONTACT  POINT 
RG-RADIUS  OF  GEAR  CONTACT  POINT 

FIND  AUXILIARY  VALUES  FOR  CALCULATION 
RP=FL0AT(Nl)/2./PN 
RG=RP’!RMPC 
C= (RP^RG) "COE 
CA-DCOS (ALPHA) 

SA-DSIN (ALPHA) 

CK-DCOS (RKS I ) 

SK-DS IN (RKS I ) 

SAK-DS IN (ALPHA- RKS I) 

C.AK-DCOS  (ALPHA-RKS I ) 

RL-RC/SA 

Al-RL-HDC/CK 

NC0EF-360.  DO-DR/FEEI /FLOAT  (ND+0.3D0 

N=NCOEF*2* 1 

NT-N-NCOEF 

FI  1  =  360. DO/RMPG/ FLOAT ( N 1 ) 

DEFINE  MATRIX  DESCIBEING  MISALIGNMENT  AND  OUTPUT  ASSEMBLING  CONDITION 
DO  5  1=1,3 
DO  5  J-1,3 
S  ( I  ,  J  ;  =0. DO 

IF  i I . EO . J )  S ( I , J ) = 1 . DO 
5  CONTINUE 
DO  15  LL-l.NG 

GAMMA-GAMMAl "FLOAT (LL-0)  /60.D0::DR 
CG-DCOS (GAMMA) 

SG-DS IN (GAMMA) 

GAMMA- GAMMA / DR  * 6  0 . DO 
IF  (NMIS.EQ.l)  THEN- 
WRITE  (LP , 10)  COE, GAMMA 

10  FORMAT  aHl,///,lX,'C-,,F5.3,'*(RP+RG);  CROSSING  ANGLE- ' , 


I 


.  %■  */ 


't 


+  F5 . 1 ,  1  (M)  ' ) 

S(1,1)-CG 
S(l,3)— SG 
S(3,1)-SG 
S (3 , 3) “CG 
ELSE 

WRITE  (LP.20)  COE, GAMMA 

20  FORMAT  (1H1.///.1X, ' C= ' ,F4.2, ' *  (RP+RC) ;  INTERSECTING  ANGLE- 
+  F5 . 1 ,  1  (M)  ' ) 

S  (2 , 2)  =CG 
S(2,3)=-SG 
S  (3 , 2)  =SG 
S (3 , 3) =CG 
END  IF 

C  SIMULATING  MESHING  OF  PINION  AND  GEAR (L=l  FOR  FINDING  INITIAL 
C  ROTATION  CORRESPONDING  TO  0  PINION  ROTATION)  AND  OUTPUT  RESULTS 
DO  25  L-1,2 
DO  35  1=1,4 
35  X(I)=0.D0 
DO  45  1=1, N 

X(7)-FEEI*FLOAT(I-(N~l)/2) 

IF  (L.E0.1)  X  C7)»0.D0 
X(5)=DARSIN(RP*X(7)*SK/R) 

X (2) =X (7) /RMPG 
SF=DSIN (X  (7) ) 

CF=DCOS (X  (7) ) 

CALL  NONLIN 
X  (8)  =X  (2)  *X  (3) 

IF  (L.EQ.l)  THEN 
XIN=X (8) 

WRITE  (LP , 30) 

30  FORMAT  (////8X,  ' FEE1 (D) ' , 8X, ' FEE2 (D) ' , 8X , ' K-ERROR  (S)  '  ,  5X, 

'RP' , 13X, 'RG'/) 

GO  TO  25 
END  IF 

X (8) =X (8) -XIN 
X  (7)  =X  (7)  /DR 
X (83 =X(8) /DR 

X  (9)  =  (X  (8)  -X  (7)  /RMPG)  *3600.  DO 
Z  (I)  =X  (8) 

ERR (I) =X(9) 

WRITE  (LP , 40)  (X(J) , J=7, 11) 

40  FORMAT  ( IX, 5F15 . 7 , F15 . 7) 

45  CONTINUE 

WRITE  (LP ,50) 

50  FORMAT  (//,'  FIND  THE  WORKING  RANGE  FOR  ONE  TOOTH:'/) 

DO  55  1=1, NT 

X(7)=FEEI;'TLOAT(I-(N*l)/2)/DR/2.DO 
X (8) =X (7) -FI  I 
KK= I^NCOEF 
ANGLE=Z(KK)-Z(I) 

ERROR ( I) = (ANGLE-FII) *3600. DO 
WRITE  (LP ,60)  X(7) ,X(8) .ANGLE 

60  FORMAT  (IX, ' (' ,F7.2, ' - ' . F7. 2, ' ) : * , F15 . 7, F15. 7) 

55  CONTINUE 


DO  65  1=1. NT 
ATEMP2=ERROR (I) 
IF  (I.NE.l)  THEN 
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IF  (ATEMP1.GT.0.D0.AND.ATEMP2.LE.0.D0)  GOTO  75 
END  IF 

ATEMP1=ATEMP2 
65  CONTINUE 

WRITE  (LP.70) 

70  FORMAT  (//IX, 'MESHING  IS  DISCONTINOUS ’ ) 

GOTO  25 

75  IF  (DABS (ATEMP1) .LT. DABS (ATEMP2) )  1=1-1 
EMAX=0 . DO 
EMIN=0 . DO 
DO  85  J=1 ,NCOEF 
KS=I+J-1 
ET=ERR (KS) 

IF  (ET.LT.EMIN)  EMIN=ET 
IF  (ET.GT.EMAX)  EMAX=ET 
85  CONTINUE 

ET-EMAX-EMIN 

KK-I+NCOEF 

WRITE  (LP.80)  Z (I) . Z (KK) , ET 

80  FORMAT  (//IX, 'WORKING  RANGE  FOR  THE  GEAR  TOOTH:  ',F7.2.' - ', 

F7 . 2/ IX, ' THE  MAXIMUM  KINEMATIC  ERROR:  ',F15.7,'  (S)‘) 

25  CONTINUE 
15  CONTINUE 

WRITE  (LP.90) 

90  FORMAT ( 1 H 1 ) 

STOP 

END 

C 

SUBROUTINE  FUNC 
C 

C  THIS  SUBROUTINE  IS  USED  TO  GENERATE  FIVE  FUNCTIONS,  THAT  IS,  AT 
C  CONTACT  POINT,  POSITION  VECTOR  AND  NORMAL  OF  PINION  AND  GEAR 
C  MUST  COINCIDE 

IMPLICIT  REAL*8(A-H,0~Z) 

COMMON  /BLOCA 1 /  X  ( 1 1 )  ,  Y  ( 1 0)  ,  A  (10 , 10)  ,  Y 1  ( 10)  , IPVT(IO) .WORK (10) , 
EPS  I , DELTA, NC , NE , NDIM 

COMMON  /BLOCA2/  ALPHA, RKSI , RP , RG , RMPG, RL , SA, CA, DR , CK, SK , HDC , 

C  ,  S  (3 , 3)  , SF , CF , SAK, CAK, Al . XP . YP , ZP , XG , YG . R 

CT=DCOS (X (4) ) 

ST=DSIN (X  (4) ) 

CFFE=DCOS (X ( 3) ) 

SFEE=DS IN (X (3) ) 

CFEK=DCOS (X (3) +RKSI) 

SFEK=DSIN(X(3)+RKSI) 

CAB=DCOS(ALPHA^X(5)) 

SAB=DSIN(ALPHA+X(5)) 

CAB2=DCOS (aLPHA-X(5) /2.D0) 

SAB2-DSIN(ALPHA*X(5) /2 . DO) 

SB’DSIN (X (5) ) 

SB2=DSIN (X (5) /2. DO) 

A2* Al  " SA!f  (CT~  1 ) 

X  ( 6 )  =  <  a2  "SaB_R"  Sb'  C  i .)  /  ( SAd  CaK-C  i  ■•Cad  'SAK,;  /  RP 
CFP=DCOS (X (67 ) 

SFP=DSIN (X (6) ) 

CAKF=DC0S(ALPHA-RKSI-X(6) ) 

SAKF=DSIN(ALPHA-RKSI-X(6> ) 

a3=2  .  D0*R’;SB2 

T3  =  CAB2"SAKF-CT’:SAB2,:CAKF 


90 


XP-A3* (CAB2*SAKF-CT*SAB2*CAKF) +A2*CAKF-RP*X (6) *CFP+RP*SFP 

YP-A3*  (CAB2*CAKF+CT*SAB2*SAKF)  -A2”'SAKF+RP*X  (6)  *SFP+RP*CFP 

ZP- (Al*SA-A3*SAB2) *ST 

RNXP=CT*CAB*CAKF+SAB*SAKF 

RNYP— CT*CAB*SAKF+SAB*CAKF 

RNZP=ST*CAB 

XG=RG*SFEE+RG*X (2) *CK*CFEK 
YG“-RG*CFEE+RG’’CX  (2)  *CK*SFEK 
ZG=X  (1) 

RNXG=CFEK 
RNYG=SFEK 
RNZG=0 . DO 
T 1  =CF*XP+SF*YP 

T2-S (1 , 1) *XC+S (1,2) *YG+S  (1 ,3)  *ZG 
WRITE  (6,50)  T1 , T2 , XP , YP , X (6) 

50  FORMAT  (IX, ’ T 1= 1 ,F15.7,5X, ‘ T2= ’ , F15 . 7 , 2F15 . 7) 

Y ( 1) =CF*XP+SF*YP-S ( 1 , 1) *XG-S (1 , 2) *YG-S (1 . 3) *ZG 
Y (2)  =CF*YP-SF*XP-S (2 , 1) *XG~S (2 , 2) *YG-S (2 , 3) *Z G-C 

Y  (3)  =ZP-S  (3  ,  1)  *XG-S  (3 , 2)  *YG-S  (3 , 3)  *ZC 

Y  (A) =CF*RNYP-SF*RNXP-S (2 , 1) *RNXG-S (2 , 2) *RNYG-S (2 , 3) *RNZG 

Y  (5)  -RNZP-S  (3  ,  1)  *RNXG-S  (3 , 2)  *RNYG-S  (3, 3)  *RNZG 
X ( 10) =DSQRT (XP*XP* YP*YP) 

X  ( 1 1 )  =DSQRT  (XG*XG+  YG*YG) 

WRITE  (6,20)  XP, YP.ZP.XG, YG,X(6) ,Al,A2,A3 
20  FORMAT  (IX, '$$$$' .6F15.7) 

RETURN 

END 

SUBROUTINE  NONLIN 

THIS  SUBROUTINE  IS  USED  TO  SOLVE  NONLINEAR  EQUATIONS  BY  NEWTON- 
RAPHSON  METHOD 
IMPLICIT  REAL’-8(A-H,0-Z) 

COMMON  /BLOCAl /  X  (1 1)  , Y (10) , A (10, 10)  , Y1 ( 10)  , IPVT (10)  . WORK ( 10) , 
EPS  I , DELTA, NC,NE,NDIM 

DO  5  1=1, NC 
CALL  FUNC 

WRITE  (6,10)  I,  (X(J) ,Y(J) ,J=1,5) 

10  FORMAT (IX, '***' , 15/5 (IX, 2D15. 7/) ) 

DO  15  J= 1 , NE 

IF  (DABS (Y(J)).GT.EPSI)  GO  TO  25 
15  CONTINUE 
GO  TO  105 
25  DO  35  J-l.NE 
35  Y1(J)=Y(J) 

DO  A5  J= 1 , NE 
DIFF=DABS  (X  ( J)  )  ’’DELTA 
IF  (X(J) .EQ.O.DO)  DIFF=DELTA 
XMAM=X ( J) 

X(J)=X(J)-DIFF 
CALL  FUNC 
X ( J )  =  XMAM 
DO  55  l\=  1 ,  NE 

55  A(K,J)=(Y1(K)-Y(K))/DIFF 
A5  CONTINUE 

DO  65  J= 1 , NE 
65  Y(J)=-Y1(J) 
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CALL  DECOMP  (NDIM, NE, A, COND, IPVT .WORK) 
CALL  SOLVE  (NDIM.NE, A, Y, IPVT) 

DO  75  J-l.NE 
75  X  ( J)  “X  ( J)  +Y  ( J) 

5  CONTINUE 
105  RETURN 
END 
C 


SUBROUTINE  DECOMP  (NDIM, N , A, COND, IPVT , WORK) 

C 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  A (NDIM, N) ,WORK(N) , IPVT(N) 

C 

C  DECOMPOSES  AREAL  MATRIX  BY  GAUSSIAN  ELIMINATION, 

C  AND  ESTIMATES  THE  CONDITION  OF  THE  MATRIX. 

C 

C  -COMPUTER  METHODS  FOR  MATHEMATICAL  COMPUTATIONS-,  BY  G.  E.  FORSYTHE, 
C  M.  A.  MALCOLM,  AND  C.  B.  MOLER  (PRENTICE-HALL,  1977) 

C 

C  USE  SUBROUTINE  SOLVE  TO  COMPUTE  SOLUTIONS  TO  LINEAR  SYSTEM. 

C 

C  IN1UT. . 

C 

C  NDIM  =  DECLARED  ROW  DIMENSION  OF  THE  ARRAY  CONTAINING  A 

C  N  ORDER  OF  THE  MATRIX 

C  A  MATRIX  TO  BE  TR I ANGULAR IZED 

C 

C  OUTPUT.. 

C 

C  A  CONTAINS  AN  UPPER  TRIANGULAR  MATRIX  U  AND  A  PREMUTED 
C  VERSION  OF  A  LOWER  TRIANGULAR  MATRIX  I-L  SO  THAT 

C  (PERMUTATION  MATRIX) *A=L*U 

C 

C  COND  =  AN  ESTIMATE  OF  THE  CONDITION  OF  A. 

C  FOR  THE  LINEAR  SYSTEM  A*X  -  B  ,  CHANGES  IN  A  AND  B 

C  MAY  CAUSE  CHANGES  COND  TIMES  AS  LARGE  IN  X. 

C  IF  COND* 1 . 0  .EQ.  COND  ,  A  IS  SINGULAR  TO  WORKING 

C  PRECISION.  COND  IS  SET  TO  1.0D+32  IF  EXACT 

C  SINGULARITY  IS  DETECTED. 

C 

C  IPVT  =  THE  PIVOT  VECTOR 

C  IPVT  (K)  =  THE  INDEX  OF  THE  K-TH  PIVOT  ROW 

C  IPVT (N)  =  (-1)** (NUMBER  OF  INTERCHANGES) 

C 

C  WORK  SPACE..  THE  VECTOR  WORK  MUST  BE  DECLARED  .AND  INCLUDED 
C  IN  THE  CALL.  ITS  IN1UT  CONTENTS  ARE  IGNORED. 

C  ITS  OUTPUT  CONTENTS  ARE  USUALLY  UNIMPORTANT. 

C 

C  THE  DETERMINANT  OF  A  CAN  BE  OBTAINED  ON  OUTPUT  BY 
C  DET(A)  =  IPVT (N)  *  A(l,i)  A  A(2,2)  *  ...  *  A(N.N)  . 

C 

IPVT (N) “ 1 

IF  (N.EQ. 1)  GO  TO  150 
NM1=N-1 

C  COMPUTE  THE  1-NORM  OF  A  . 

AN0RM=0 . DO 
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DO  20  J-l.N 
T=O.DO 
DO  10  1  =  1, N 
10  T=T+DABS(A(I, J)) 

IF  (T.GT.ANORM)  ANORM-T 
20  CONTINUE 


DO  GAUSSIAN  ELIMINATION  WITH  PARTIAL 
PIVOTING. 


DO  70  K= 1 , NM1 
KP 1=K+ 1 


FIND  THE  PIVOT. 


30 


M-K 

DO  30  I=KP1,N 

IF  (DABS (A(I ,K) ) .GT.DABS (a(M,K) ) )  M=I 
CONTINUE 
IP VT  (K) =M 

IF  (M.NE.K)  IPVT(N)=-IPVT(N) 

T-A(M.K) 


A  (M .  K)  =A  (K ,  K) 
A(K,K)-T 


SAIP  THE  ELIMINATION  STEP  IF  PIVOT  IS  ZERO. 
IF  (T.EQ.O.DO)  GO  TO  70 


COMPUTE  THE  MULTIPLIERS. 


AO 


DO  AO  I=KP 1 , N 
A(I,K)=-A(I,K)/T 


INTERCHANGE  AND  ELIMINATE  BY  COLUMNS. 


50 

60 


DO  60  J-KPl.N 
T=A(M. J) 

A (M, J) =A (K, J) 

A(K,J)=T 

IF  (T.EQ.O.DO)  GO  TO  60 
DO  50  I=KP 1 , N 
A(I, J)»A(I,J)+A(I,K)*T 
CONTINUE 


70  CONTINUE 


COND  =  (1-NORM  OF  A) * (AN  ESTIMATE  OF  THE  1-NORM  OF  A-INVERSE) 

THE  ESTIMATE  IS  OBTAINED  BY  ONE  STEP  OF  INVERSE  ITERATION  FOR  THE 
SMALL  SINGULAR  VECTOR.  THIS  INVOLVES  SOLVING  TWO  SYSTEMS 
OF  EQUATIONS,  (A- TRANSPOSE)  *Y  =  E  AND  A::Z  =  Y  WHERE  E 

IS  A  VECTOR  OF  +1  OR  -1  COMPONENTS  CHOSEN  TO  CAUSS  GROWTH  IN  Y. 
ESTIMATE  =  (1-NORM  OF  Z)/ (1-NORM  OF  Y) 


SOLVE  (A-TRANSPOSE) *Y  =  E 


,  N 


DO  100  K= 1 , 

T=0 . DO 

IF  fK.EQ. 1)  GO  TO  90 


80 

90 


KMl'K-1 
DO  80  1=1, KMl 
T=T-A ( I , K) "WORK (I) 

EK=1 .DO 
IF  (T.LT.O.DO)  EK=- 1 . DO 
IF  (A(K,K) .EQ.O.DO)  GO  TO  160 
100  WORK(K) =- (EK+T) /A(K,K) 

DO  120  KB=1 ,NMi 
K=N-KB 
T=0 .  DO 
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M 
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KP1-K+1 

DO  110  I-KPl.N 
110  T=T+A(I,K)*WORK(K) 

WORK (K) “T 
M=IPVT(K) 

IF  (M.EQ.K)  GO  TO  120 
T-WORK(M) 

WORK (M) “WORK (K) 

WORK (K) =T 
120  CONTINUE 
C 

YNORM'O . DO 
DO  130  1=1, N 

130  YNCF>M=  YNOPM+DABS  (WORK  ( I)  ) 

C 

C  SOLVE  A*Z  =  Y 

CALL  SOLVE  (NDIM , N , A , WORK, IPVT) 

C 

ZNORM=0 . DO 
DO  140  1  =  1, N 

140  ZNORM=ZNORM+DABS(WORK(l)) 

C 

C  ESTIMATE  THE  CONDITION. 

COND=ANORM,fZNORM/ YNORM 
IF  (CON'D.  LT.  1 .  DO)  COND=l.D0 
RETURN 

C  1-BY-l  CASE. . 

150  COND-l.DO 

IF  (A (1 , 1) .NE.O.DO)  RETURN 
C 

C  EXACT  SINGULARITY 

160  COND= 1 . 0D32 
RETURN- 
END 

SUBROUTINE  SOLVE  (NDIM, N , A, B , IPVT) 

C 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  A (NDIM, N) ,B (N) , IPVT(N) 

C 

C  SOLVES  A  LINEAR  SYSTEM,  A*X  =  B 

C  DO  NOT  SOLVE  THE  SYSTEM  IF  DECOMP  HAS  DETECTED  SINGULARITY. 

C 

C  -COMPUTER  METHODS  FOR  MATHEMATICAL  COMPUTATIONS-,  BY  G.  E.  rORSYTHE , 
C  M.  A.  MALCOLM,  AND  C.  B.  MOLER  (PRENTICE-HALL,  1977) 

C 

C  INIL'T .  . 

C 

C  NDIM  =  DECLARED  ROW  DIMENSION  OF  ARRAY  CONTAINING  A 

C  N  ORDER  OF  MATRIX 

C  A  TRIANGULAR  I ZED  MATRIX  OBTAINED  FROM  SUBROUTINE  DECOMP 

C  B  RIGHT  HAND  SIDE  VECTOR 

C  IPVT  =  PIVOT  VECTOR  OBTAINED  FROM  DECOMP 

C 

C  OUTPUT.. 

C 

C  B  =  SOLUTION  VECTOR,  X 
C 

C  DO  THE  FORWARD  ELIMINATION. 
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IF  (N.EQ.l)  GO  TO  50 
NM1=N-1 
DO  20  K-l.NMl 
KP1=K+1 
M=IPVT(K) 

T=B  (M) 

B(M)=B(K) 

B  (K)  “T 

DO  10  I-KPl.N 
10  B(I)=B(I)+A(I,K)*T 

20  CONTINUE 

DO  40  KB= 1 , NM1 
KM1=N-KB 
K=KM1+1 

b(k)=b(k)/a(k,k) 

T— B(K) 

DO  30  1=1, KM1 
30  B(I)-B(I)+A(I,K)*T 

40  CONTINUE 
50  B(1)=B(1)/A(1,1) 
RETURN 
END 


NOW  DO  THE  BACA  SUBSTITUTION 


iV.  A-V.  A  a  -V.  A  •  , 
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